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A MICROCOMPUTER-BASED POSITION UPDATING SYSTEM FOR 
GENERAL AVIATION UTILIZING LORAN-C 


This report Identifies the advantages Loran-C offers over 
VOR/DME as a primary navigation aid. These advantages Include 
long range coverage by several stations; coverage at low alti- 
tudes and the capability to have non-precision approaches at 
airports not already served by landing aids* Modem digital 
electronic technology is used to produce a device to convert 
Loran-C to useful pilot Informarlon using simple software algo- 
rithms and low-cost microprocessor devices* The cost and lack 
of availability of suitable processors to execute these algo- 
rithms have prevented a Loran-C navigator from being developed 
prior to this work* Results indicate that the microprocessor- 
based Loran-C navigator has an accuracy of 1*0 nm or less over 
an area typically covered by a triad of Loran-C stations and 
can execute a position update in less than 0*2 seconds* The 
system has been tested in 30 hours of flight and has proved 
that it can give reliable and accurate navigation Information* 
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1. INTRODUCTION AND STATEMENT OF FROBLEM OF rw' 

Loran-C, a 100 KHz, tlaa-dif farence, hyperbolic radio navigation 
acheae, haa certain advantagea for uae aa a primary navigation aid for 
general^aviation aircraft* Significant advantages include long-range 
coverage by using a few craneaitte 's, useful signals at low altitudea and 
in mountainous terrain, and signals which are stable allowing accurate 
and consistent positional data to be obtained. These advantages have 
prompted federal authorities to consider this system as a primary naviga- 
tion aid and for application as a non-precision (two-dimensional) landing 
aid. I 


Lcran-A, which ope.'atea at 2.5 MHz, was in use for over 30 years 
but has now been discontinued. Loran-C, which has been in operation for 
20 years, is continuing to expand with new chains in all parts of the 
world. Loran-C provides more stable signals than Loran-A over long 
distances and for this reason, there has been increased interest in land- 
based use of Loran-C for aviation activities. Loran has traditionally 
been used by the maritime services and has not been accepted previously 
in aviation circles because of users' unfamillarity with hyperbolic 
systems and their impracticality and the lack of low-cost, simple equip- 
ment to convert the hyperbolic data to a more readily understood form. 

This paper documents a rather complete effort to design, build, 
and flight test a navigator unit which converts Loran-C data at the out- 
put of a receiver to latituce/longitudc coordinates and also range/ 
bearing coordinates to a waypoint. This work has taken advantage of 
contemporary microprocessor technology for simple implementation of new 
algorithms to give accurate positional information from Loran— C. This 
has ii..olved investigation and testing of various data conversion 
methods, designing and building the necessary microcomputer hardware and 
the writing of the data conversion software. Finally, the entire concept 
was tested during approximately 30 hours of flight in a Piper Cherokee 
end a Douglas DC-3. It was proved that a simple method will Indeed con- 
vert Loran-C data to a pilot-usable fora. The microcomputer equipment 
necessary is low in cost and readily available and the system works %«ell 
in providing good navigation information in an aircraft flight environ- 
ment. With the continuing trend toward more powerful and lower cost 
electronic devices, the roadblocks to implementing fully a Loran-C navi- 
gator have essentially been removed. 

The international standard csethod of navigating is to use the VOR 
(VHF Omni Range) snd DME (Distance Measuring Equipment). This system, 
with over 1000 stations in the United States, has proven to be acceptable 


^Non-precision is used to connote that no vertical guidance is pro- 
vided. An extensive study in this area has teen ongoing in the 
state of Vermont. See, for example, Polhemus, W.L. , ”Evaluatlon of 
Loran-C for Enroute Navigation and non-Precision Approach Within 
the State of Vermont," Proceedings of the National Aerospace 
Symposium, Dayton, Ohio, March 1980. 
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and very capable. InforoaClon la provided to the pilot In the forw t.f 
the angular bearing and distance relative to a VOK/DME station [1]. A 
more recent feature called area navigation (RNAV) works in conjunction 
with VOR/DME d.:ta to allow the pilot to set up arbitrary waypoints (at 
airports, VOK/DHE stations, etc.) and then fly a stralght-llne course 
between these points. 

The incentive to consider low-frequency systems is that the VHF 
and UHF signals from a VOK/DME station do not propagate to follow the 
earth and provide low altitude coverage and signals and can be blocked by 
mountains. Also, with low-frequency systems, 20 to 30 stations could 
cover the entire continental U.S. which would be considerably cheaper and 
easier to maintain than over 1000 stations. 

With the ability to cover low altitudes, Loran-C could provide, 
at no extra cost, an accurate approach navigation system currently un- 
available at many airports. There are over 14,000 airports In the U.S., 
half of %ihlch do not hive electronic instrument landing aids. By having 
an all-weather approach :ld, such as Loran-C could provide, the utility 
and cost effectiveness of existing airports could be Increased. 

Clearly, there are other methods of air navigation Including non- 
dlrectlonal beacons (NDB), dead reckoning or inertial navigation. There 
Is also some promise of using satellites for navigation [2]. These 
methods have disadvantages. For example, the NDB is not as precise 
to use as the Loran-C and Inertial navigation and satellite systems are 
certainly much more expensive than Loran-C. 

The present-day, Loran-C network currently covers much of the 
northern hemisphere as shown In Figure 1-1. The possibility of using 
this system especially for non-precision approaches and RNAV has 
generated much Interest and has been investigated extensively [3,4]. 

This research has dealt with ways to handle such practical problems as 
GOOF (geometric dilution of precision), which occurs when the hyperbolic 
llnes-of-posltlon cross at acute angles, and the effect of passing over 
a Loran-C station. It has been shown that these problems can be handled 
effectively without seriously affecting the overall utility of the system 
[5,61. 
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This brief review of low»frequency navigation «y a terns, 

Including Loran, la intended to show how present-day Loran-C has bene- 
fitted froa experiences gained through experlaents with ocher syutums. 
The operation of the Loran-C systea is discussed, and the details of 
how navigation inforaation is derived froa Loran-C are preaented. 


A. Historical Developaent . Early in aviation history it 
bccaae apparent Chat flying activities needed to be conducted at claen 
when navigation by visual cues alone was not adequate. Other iseans oC 
navigation weie needed, and these can be divided into self-contained 
and ground-based systeas. Self-contained systeas such as Inertial 
navigation have all the navigation conponents on board the aircraft. 

A ground-based systea sends a signal or signals to Che aircraft where 
navigation Inforaation is derived. Ground-based systeas, such as 
VORs, are preferred because the heavy and bulky equipment can be kept 
off Che aircraft. In addition, use of ground-based ntvlgatlon systeas 
can provide correlation among all users to a fixed reference point 
(or a set of points fixed with respect to the users) thus allowing 
users to establish their position and avoid collisions with each 
other. 


Reasons for Low-Frequency Systeas . It ir aost "^tural 
to use radio slgnala to transalt navigation inforaation to an 
aircraft, although there are restrictions as Co which frequencies are 
opclaua. As shown In Table 2-1 the propagation characteristics 
of various frequency regions vary considerably [7,8]; at very low fre- 
quencies the transaltted waveform follows the earth's curvature In a 
duct foraed between the earth and the Ionosphere, and in the VHP 
region the waveform follows a straight-line path from the transmitter 
to the receiver and can be blocked by the horizon or large objects 
such as iBountains. 

While low-frequency signals may be transmitted over several 
hundreds or thousands of kiloaeters, with sufficient power, these 
signals often are contaainated with noise to the extent that the level 
of average noise Is higher Chan the signal. The transmission of 
aedlua-frequency signals for navigation use has a serious disadvantage 
in Chat these signals are refracted by the ionosphere causing the path 
length of Che signal to cK^nge as the height of the Ionosphere v'arles 
with clae of day ^nd season. The use of VUF or higher frequencies 
affords btndwldchs of several kilohertz to Mgahertz lor transmitting 
navigation inforaation; however, with llne-of-sight coverage, the cost 
of placing and maintaining the transmitCers necessary ‘.o provide adequate 
coverage la a serious disadvantage. 

For covering a large area of several hundred thousand kilome- 
ters or aore, a low-frequency systea appears appropriate. At very low 
frequencies, transmitter and antenna design Is difficult because of 
the high power levels and long wavelengths involved. The bandwidth 
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at vary low fraquanclaa la narrow, laas than Chat rwqulreJ 
for voice cranaaiaaton, thua Halting ^.ow auch Information may t>e 
aant [9] . The low-fraquancy band allowa enough bandwidth to u.ie 
pulaa-’Cyp^' navigation ayataaa which alapllfy '.atection of navigation 
Information coaparad to the uae of con Clnuour -wave ayatems. 

2. Con> rlbutlona of Other S yatema . Many CeaCa were con- 
ducted ualng low ffaquenclaa ^or navigation. Some of the teata led Co 
ayateaa that did worh., although they were not very practical and were 
difficult to uoo. Becauae of the stable nature of very-Iow-fiequency 
signals, some systems were developed which used variations In pliaac 
ralatlva to receiver location. The use of continuously transalcted 
signals was nacesaary because of the bandwidth restriction on the 
antennas. The European Dacca and the more recently developed Omega 
ays teas are exaaples of these. In these systeas, the receiver com- 
pares Che relative phase between two or store stations, ar.J these phase 
differences Chen define a set of hyperbolic lines of position, the 
Intersection of which Is the receiver's location. A lane Iden- 
tification probleEi Is Inherent to these systems because of the cyclic 
redundancy of the transalcted signal. There can be considerable 
variation of Che transalcted phase because of natural phenomena, and 
this has been the subject of such research [10,11]. 

In the 1940s, work was begun on pulse-type navigation systems 
operating at various frequencies, mostly In the mcdlum-’f requency 
range [12]. These systeas attempted Co derive navigation Information 
by uilng Che time difference of arrival between pulses received from 
different scat lone. Over a period of tlswe, the technology of 
designing highly stable oscillators became very refined and It was 
soon recognised chtat siuch better navigation coverage could be obtained 
In Che low-frequency band. Because of the demands of Che Second World 
War, a pulse-type system, called Lo.ran-A and operating near 2 MHs, was 
established while research on low-f requsney techniques continued. In 
1946, a system callec Cyclau using automatic cycle Identification and 
phase measuring techniques was developed. This soon led to the 
testing of another system called Cytac, operating In the 90-100 kHz 
band, which used the leading edge of pulses to neasure time differ- 
ences. In the late 1930s, a set of specif lea t'.ons was developed which 
was considerably improved over Che standard Loran, and thus the Cytac 
tests became what Is known today as Loran-C [13]. The Loran-A system 
was Chen phased out beginning In the late 1970s with complete 
withdrawal of the system In 1980.^ 

Throughout many early experiments there was very little known 
about Che propagation characteristics at low frequencies and how these 
characteristics would Influence Che outcome of the navigation experi- 
ment. Often, some rather surprising discoveries were made during the 


^A good discussion of these riavlgation methods is provided by 
Sandrettu, P.C., Aeronautical Radio Engineering , McCraw-hlll Book 
Company, New York, 1942. 
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testing of these prospective systems and iiach knowledge was gained by 
actual field «aork. Loran-C today contains many elemoncs such as long* 
range covers, ^e and pulse transmission to avoid ositual interference and 
skywave conttmii^at on, which were discovered years earlier. The 
system is mr. optimization of many of its predecessors. 

P. Basics of Loran’-C Time Differences (TD) . Each Loran-C sta- 
tior transmits a serieo'of pulses with a carrier frequency of 100 kHz. 
An examplv. of a transmitted pulse is shown in Figure 2-1. The r. the- 
matical description of this pulse is: 


f(t) -( I" r.xp[2(t - -p) J Cosu^t 
P P 


%ihcre tp “ 65 and - 2vxI00 kHz. Thi i equation repreaenta a 
pulse with rapid rise tine to allow the receiver to lock onto the 
early part of the pulse before a sk\wave pulse appears. The latter 
part of the pulse has a slow decay time to confine the energy to a 
band extending from 80 to 120 kHz. fhe phase of the leading edge Is 
tightly controlled since this is where the receiver (siamples the signal 
to determine time differences [14]. 

I. Development of Time Dlffereaces . Loran-'C stations 
associated in groups of ^>s tween three ancT five stations; one station 
is designated the master and the others are designated secondary or 
slave stations. As shown In i*lgare 2-2. the master station transmits 
nine pulses, the first eight separated from each other by 1000 ps, and 
the ninth separated from the eighth by 2000 ps. This Is done for 
identification purposes. Each of the secondary stations transmits 
eight prilses separated from each other by 1000 ps. A particular i:roup 
of station* is called a chain and is identified by the chained group 
repetition Interval (CRl) which is the rate that a station transmits 
pulse groups. Further identification is provided by periodically 
reversing the phase of the carrier for each station in a set sequence 
(151. Typical GRIs are in the range of 49900 ps to 99900 ms. By 
measuring the amount of time between the arrival of a pulse set for 
different stations, the Loran-C time differences are determined which 
can then be used to 1< cate the receiver's position. This will oe 
examined in more detail in the next chapter* 

2 Hype rb olic Coordinates . As will be pc inted out later, 
the Loran-C TD readings are of the form of a hyperbolic equation. 

This places Loran--C with most other navigation systems used below 
30 HHx whose primary source of information is a set of hyperbolic 
coordinates. Ihe reasou for this Is that narrow bandwldths at low 
frequencies do not allow sufficient data rates to transmit slgnelc 
N>lth a higher information content than the pulse system. Unfortu- 
nately, it also puts more of a burden on the receiver or operator who 
must then convert the coordinates to a mere convenient form. The most 
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Loran-C Trc'-smlfred V/avefor.r.. Frequency of corrler is 100 KHz. 





natural form of navigation information is the user's location is rela- 
tive to a reference point* A hyperbolic system only provides an 
abstract set of coordinates defined by the particular setup of sta- 
tions* 


3* Llnes^f-Fosltlon * By measuring the difference in time 
of arrival of a set of pulses from two stations* a TD number is found. 

Referring to Figure 2**3» a TD number is a point on the locus cf points 

all having the same value. Denoting the distances from each of sta- 
tions to the point as d. and d 2 * and the constant difference as 2a* 

then 


and 




-h y^ 
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(2-1) 




x-c) + y 


( 2 - 2 ) 


Since the difference of the distance is constant* i.e.* ^2 ” 

then 


/ 


(x4c) + y 


-C. 


y)2 + 


2a 


(2-3) 


Rearranging and squaring both sides 


(x4c)^ + + 2a + y^ + (x-c)^ + y^ (2-4) 

By further rearranging equation (2-4) 

(a^ - c^)x2 - a^y^ - a^(a^ - c^) (2-5) 

2 2 2 

substlcutlng c ■ a + b . the standard hyperbolic equation is pro- 
duced. That is ' ^ 



1 


( 2 - 6 ) 
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Point a r<^*|>r*-s»!ucs the x-axis crossing of the hyperbola and 
point b is Che conjugate axis l»^ngCh such that the asymptote of the 
hyperbola is y - 

a 

Each hyi>erl>oirc locus of points is a llne-of-posit ion (LOP). 
Since one LOP does not uniquely define the receiver's position, a 
second master*secondary pair roust l>e used to define a second l^P; the 
intersection of two or laore LOPs defines the n!ceLver*s position, as 
shown in Figure 2-4. 

^tethod of computing time-di fferences. Reference for the 
following discussion will be made to Figure 2-5. Tne signal 
transmitted by the master station (M) is radiated in all direct I ons. 
The time it takes for the signal to reach a given point is related to 
the S|H 2 ed of light (approximately 300 meters per ys) and to a lesser 
extent, the medium through which the signal travels. The secondary 
stations always transmit after the roaster; the time interval between 
the transmission of the master signal and the secondary signal is 
equal to the amount of time it would take the master's signal to reach 
the secondary (baseline time) plus a coding delay. Coding delays are 
Introduced at the secondaries to ensure that no two stations will 
overlap anywhere in the coverage region. At the receiver position, 
there will be two signals present, the roaster and the secondary 
separated by a certain amount of time. The amount of time can be 
expressed ma thematically as: 


TD - B + A + Tg - Tj^ 


(2-7) 


where TD is the time difference in arrival of the master and secondary 
signals, B is the one-way baseline time between the troister and secon- 
dary, A is the secondary coding delay inserted by the secondary, T 3 is 
the one-way baseline time between the secondary station and the 
receiver, and T)^ is the one-way baseline time between the receiver and 
the master station. The baseline time T|yj is subtracted from the above 
equation to account for the receiver's displacement from the master 
station. Equation (2-7) contains the difference of x and y terms in 
the form of Tg and T^^ plus additional constants and is thus in the 
form of a hyperbolic equation. 

The baseline quantity B and the coding delay are generally 
known from the values set up during the installation of the Loran-C 
chain. The baseline quantity B can be accurately measured and the 
coding delay is chosen to provide reasonable distribution of the 
secondary signals throughout the group repetition interval. These two 
quantities are usually combined and published in the data for a i^ar- 
ticular chain by the U. S. Coast Guard [16]. The two baseline times 
from Lhe receiver, Tg and are unknown and must be calculated 
before they can be applied to the TD equation. This Involves computing 
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Figur« 2*4. TOs iteceiv«c> from fh« Voriot’i MostafSecondor/ Pain Define 
LOf^ which All Intenect the Receiver's Position. Coch 
LOP corresponds to a porticulor TD nwmoer (in microseconds). 
The LOPs are actually hyperbolas with their corresponding 
moster-secondory pair of statiortt as foci. 


- 13 - 


OF 


^ /s 


'’oos Qmuty 


Secondary 



Receiver M Master 


Fi0wfe 2-5. The TO KeodinB ot the Keceivur is Formed by the Amount of Time 
the Signol TroveU Between the Mooter and SeconJory, Plu^ the 
Coding Delay, Plus the Amount of Time the Signoi TrovoU from 
the Secondory to the Receiver, Minus the Amount of Time the 
Signoi Traveled from the Mostor to the Receiver* 


- 14 - 


the length between the receiver and the Loran-C stations and computing 
how much time it would take for the Loran-C signal to traverse this 
distance* Being able to do this is the most important part of a TD 
prediction algorithm* 

Computing the amount of time required for an electromagnetic 
wave to travel from one point to another involves two factors: calcu- 
lating the exact baseline distance between the two points taking into 
account: the non-*spherical nature of the earthy and applying correc- 
tions to the velocity of the signal as it passes through the medium 
around the baseline* The conductivity of the earth and effect of the 
air have the greatest impact on the velocity of propagation. 

A method of calculating baseline length between two i>oints on 
the surface of the earth will now be presented* As sho%m in Figure 

the earth is approximated by an ellipsoid with major (equatorial) 
radius* a* and minor (polar) radius* b* The usual coordinate frame 
for defining a point on the surface of the earth ic the geodetic 
spherical coordinate system where the latitude is represented by ^ and 
longitude is represented by X. The geodesic arc-length between two 
points on the earth may be computed by the following procedure which 
calculates the arc-length on a plane through the center of the 
ellipsoid [17]. One begins by defining ^ as the geodetic latitude of 
the receiver and as the geodetic latitude of the second point 
(Loran-C transmitter) and AX as the difference in longitude between 
the receiver and second point* The procedure for computing the arc 
length is to consider a sphere of radius a circumscribing the earth. 

By projecting the latitude of a point on the earth onto the sphere* 
the point on the sphere has a parametric latitude denoted by 3 [18]. 
The parametric or reduced latitude is computed as: 


tan 3 

- (1 

~ f)tan« 

(2-8) 

can3£ 

- (1 

- f)tan^j^ 

(2-9) 


where the subscript i represents the second point and f is ^he flat- 
tening of the ellipsoid: 

f . a " b (2-10) 

a 

The generalised direction cosines of the projected point are: 


- cosS^sinAX (2-11) 

C 2 ■ cos3sin3| ~ sin3co83|Cos( AX) (2-12) 
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(2-13) 


OKlGiMHi- r> 


■■ sln3Hln0|^ •*■ cos3cos|)^co8( AA) 


The bearing angle, vlch respect Co true north, fron Che receiver to the 
second point. Measured at Che receiver is: 


tanT • — 

Co 


(2-lA) 


The approxlaaCe angle from Che receiver co Che second point. In a 
plane through the center of Che ellipsoid Is: 


CanO 


C2Cosy + Cjslnt 


(2-15) 


The geodesic arc length between the receiver and the Loran-C 
transffllCCer suty be found from 


P 



(mu + nv) 


(2-16) 


where: 


in 


(sln0 + slnpj^) 


2 


n 


sinB - sin3. 

r i ]2 

sinO ^ 


It • f ^ “ c obQ ^ ( 0 - slnQ 3 
^ slnO ^ ^ s^fnO ^ 


(2-17) 


V - (1 + cosO) (0 + slnO) 


To calculate Loran-C tine dif ferencec , it 1b necessary to do this pro- 
cedure twice; once for the receiver to master station arc length and 
once for the receiver to secondary station arc length. The units ot p 
depend on the units used to specify the radius of the earth, a. If 
the units of a are in microseconds (300 meters - I us), then p gives 
directly the amount of time needed for the Loran-C signal to travel 
from one point to another. This method uses an oblate spheroid (sodel 
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of Ch«i MrtU Co approxiiMCe Che goodoelc. Over UlttC^ittcoK typical v>t 
Cho«e involved with Lor«ri*»C uork^ Chie uaually la adequate. In a 
generalised developnent of geodesics on an ellipsoid by liolmstrom, a 
three-dleenelonal aiodel of the earth was developed allowing the three 
axes to be unequal. Using an oblate spheroid case where the equa-* 
torlal axes are equal, the denonstrated error over a 482 kis geodesic 
was 0.82 a with a 0.033 second of arc error [19] • 


Although the above procedure aiay be used to calculate Loran-C 
tlae differences, the assuaptlm made is that tM l/>ran-C signal will 
travel over the baseline path with a constant velocity, independent of 
the adjacent MdluB. In practical Loran*C «x>rk, the Loran-C signal 
travels over esdia which can have widely varying characteristics. 
Except for all-water paths, the signal travels over a surface of 
Inhottogeneous conductivity and dielectric constants and also one which 
has irregular surfaces. There have been many attempts to produce a 
working model of the earth's surface In order to resolve the problem 
of predicting the signal phase delay, but such attempts are very dii- 
ficult in light of the non-spherlcal and the irregular surface Impe- 
dance nature of the earth. In an Integral equation model of a 
nonhomogeneoua, Irregular earth proposed by Sanaddar, predicted secon- 
dary phase delays of 2.5 to 3.3 ps over a path length of 300 co 400 km 
were calculated which compare favorably to measured values [20]. The 
difficulty in applying such saodels lies In getting an accurate profile 
of the section of the earth over which the signal will |^ss. The 
numerical calculations Involved in solving the integral equations are 
complex and yield a phase factor correction of only a few microse- 
conds More general impedance models are typically «siiploycd, sucli as 
with the Defense Mapping Agency^ 


The total time difference equation including secondary phase 
factor corrections discussed above is: 


TD - ^ (^S " V ^ ^ (2-18) 

where and are the calculated baseline lengths to Che U^ran-C 
secondary and master stations, respectively, using equation (2-16) » 

V Is the index of refraction over the geodesic, C is the tree-space 
speed of light, L Is the baseline length between the ouster and secon- 
dary stations, A Is the secondary coding delay, and cs the 

secondary phase factor corrections for the secondary and master geode- 
sics, respectively. 


^The Defense Happing Agency computes corrections for overland phase 
retardation based on a data base of calibrated stations. 
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Appendix I Llatx « FQKTKAN^IV prograa run on un IHM 37C/1S8 to 
calculate Loran-C tlae diffaraucea. Varloua conatantH, Including the 
geodetic latltudaa atwl longttudca of the Loran-C aaeter and eecondary 
atatlona, the naater-to-aer.ondary baaeline tlcw delays, the secondary 
coding delays, and earth constants are all stored Internal to the 
prograa. Ihe program accepts th« coordinates of the receiver as 
Input. The oblate spherical earth model discussed above la used, 
however, no secondary phase factor corrections are applied. Figure 
2-7 shows the results of using this program with several points a’.ong 
with actual time difference values calculated for these points using 
over-land phase retardation corrections [21]. This shows a typical 
error of two to three microseconds. The main source of error In this 
program Is the lack of Input for the secondary phase factor correla- 
tions. It should be noted that the Defense Mapping Agency maintains 
phase factor corrections for only a few control points; thus their 
calculations are not always exact for the user's position. 
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III. CONVERTING TO* TO K)SIT10N OF POOR QUALITY 

The preceding chapter has dealt with how to cnnpute the Loran-C 
tloM* differences which My be expected at a given location* Aa will 
be seen, this is necessary for soiae iterative U>raa Tl>~to-posit Ion 
solutions* Now, the Inverse problem will he considered, that of con- 
verting received Loran-C time differences to more generalized position 
coordinates, such as, latitude and longitude. In this chapter, 
several different techniques will he examined to evaluate their 
simplicity, accuracy, and suitability for being used In a small 
microcomputer-based navigation system. Some methods of converting 
rhesc first-order derived position coordinates Into a range and 
bearing form simulating a VORTAC reference which Is more familiar to 
pilots will also be discubsed* 

A* Conversion Methods Available . The time difference readings 
obtained from the Loran-C receiver are generally of lUtlo use to a 
pilot* Even if a pilot could become familiar with the time differ- 
ences in certain areas through experience, using similar time dif- 
ferences for a different triad of stations or with a dlfl^rent chain 
would be meaningless* As shown in Figure >-l, there Is no simple 
relationship between the Loran-C llnes-of-posltion and a geocentric 
grid coordinate system, such as latitude-longitude. Further, the 
spacing between successive LOPs varies with the position ^ the cover- 
age region; at large distances between stations, a change of 10 ps 
might relate to a linear change of 10 km, while at a distance close co 

the station, a change of 10 ps might correspond to a linear change o^ 

1 Km. Thus, navigating by Loran-C time differences alone is lot 
a very suitable method* 

The classic method of relating tlirj differences to a more uni- 
versal set of coordinates is through the use of charts or cables* 

While this My suffice for maritime services (for which such charts 
and maps are Mde) where craft velocity is low, it is very unsuitable 

for aircraft use. Ideally, the navigation unit should allow t>»e pilot 

to Instruct Che unit what course the pilot wants to fly, and Che unit 
will then provide the pilot with course gu ..dance Information as the 
flight progresses without the pilot iiavlng to derive course guidance 
from M^:^ or tables* For this reason, it is desired to calculate 
autoMtically the actual posltic:n in a universal set of coordinates 
from Che Loran-C time differences. This information may then be used 
for input to a course director computer. The course direction com- 
puter would then drive a display indicating whether or not the 
aircraft is on course, and the pilot wot;i<l then simply maintain an on- 
course indication* 

There are Mny means available for converting U)ran-C time dif- 
ferences to latitude and longitude. These may be divided into two 
basic processes, those which require repetitive calculations to arrive 
at the final result, and those which us^i ^ direct or closed-form solu- 
tion to do the conversion. 
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An iterative TD-to-poblt Ion algorithm functlonB by assuming 
that the receiver la at a given location, then comparing the ixisltlon 
determined by the received TDs ti the assumed position. The TUs lor 
the assumed position are calculated before the comparison and then the 
assumed position is adjusted so as to minimize the error In comparison 
of the two sets of TDs* A way of mechanizing this U shown in Figure 
3-2, where for a given triad of Loran-C Btations, a table of TD values 
Is generated spanning the coverage region of the triad, and another 
table is generated to the corresponding positions for the TD table. 
When a set of TDs are received, a process of linear interpolation is 
done between the two tables to find the actual position. This could 
be done several times to increase the accuracy of the fix. Other ways 
of relating received TD values to assumed TD values will be discussed 
in Chapter Xli*C. 

A non-iterative, or closed-form solution, calculates the pre- 
sent position of the receiver by setting up a set of spherical 
equations In such a way tliaC the received TDs can be used as parame- 
ters in these equations, then the equations are solved for the unknown 
values representing the receiver's position. This is generally 
complicated because spherical equations have angles as their argu- 
ments, and relationships between position on a geocentric grid and 
these spherical angles must be found. Once this is done, the 
equations are solved, generally with a computer, and the resulting 
functions are related back to the geocentric grid to find the latitude 
and longitude of the receiver corresponding to the received TUs. 

B. Analysi s of Conversion Methods . There are many specific 
TD-to-posltlon methods to choose from, and there raay be several that 
perform In a similar fashion. At this point, several criteria will bo 
developed Co select one for use In the microprocessor controlled navi- 
gator unit. 

The overall accuracy of any particular method Is important; 
some conversion methods are simple from a computational point of view, 
but suffer from poor accuracy and possible instabilities in their 
calculation. For aircraft use, Che conversion method should have an 
accuracy of better chan one kilometer; accuracies on the order of 
0*1 nm should be more typical. This assumes that the received TDs 
correspond exactly to the actual position; then the TD-to-poslc ion 
conversion will converge to this point with a possible error of 
O.l nm. This confidence interval assures that a pilot will be able to 
reach a given navigation fix or airport with good accuracy. 

A second limiting factor on choosing a conversion inethod Is how 
much computational time will be required. For maritime use, where 
craft velocities are low, computation times of one minute, or more, 
are reasonable. However, for aviation purposes, a considerably faster 
position computation is needed. Host general purpose microprocessors 
on the market today, do :iot lend themselves very well to arithematlc 
calculations; typically, there are many involved calculations that 
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fiecd to bo made la a TD-tO'-poslC loa coavorbLon. Thu«> a coiwerstua 
alogrlchn should involve a minimum of compIcK type mathematical calcu- 
lations , particularly transcendental functions^ so that the micropro- 
cessor system can perform a TD-to-poal t ion conversion in less than one 
second. 


Finally, the size of the conversion algorithm needs to be con- 
sidered. In most microprocessor applications, it is desirable to keep 
the memory size low to avoid excessive power drain and other opera- 
tional problems that can arise in large memory systems. Even on large 
computer systems, a small conversion program is easier to store, and a 
program that requires a small amount of memory i easier to manage in 
main storage and will often run faster for that reason alone. 


An Algorithm for Hicrocomputer 'Jse. Several possible 
TD-'to-posltion conversion methods will now be analyzed using the con- 
cepts developed above* Iterative conversions are examined first. 

As mentioned briefly in Chapter Ili-l, after receiving a i*air 
of Loran-C time differences, it is necessary to compute the amount and 
direction in which the initial position should be adjusted. A pro- 
perty of hyperbolic navigation systems is helpful In this problem 
(Figure 3-3); referring again to the TD equation TD • 0-l-A+Tc^-Tj^, It 
can be seen that as the receiver is moved closer to the nkister 
(farther from the secondary), will increase and will decrease, 
and the TD value for that master-secondary pair will increase. Thus 
if the difference between the assumed position and the measured TDs is 
positive, the new estimate of position should be made closer to the 
master station to increase the value of the calculated T!) and decrease 
the result of the subtraction of the two TDs. This process is 
illustrated in Figure 3-4. 

Note that this algorithm could take many Iterations to reacli an 
acceptable estimate of actual position. The amount of position change 
to be made with respect to the magnitude of the TD error must be care- 
fully chosen for the chain and particular stations used. Depending on 
the actual position relative to the lx>ran-C stations and any b.isellne 
extensions between stations, this routine could suffer from a large 
radius of convergence and poor convergence time. 

Another method of calculating the actual position based on the 
received TD» is to assume an initial reference point and calculate the 
TDs at that point as before. Then the measured TDs are compared witli 
the TDs based on the assumed position. T^e assumed position is then 
moved to a new position based on the magnitude and sign of the TD 
error. This procedure differs from the previously discussed case in 
that an attempt is tsade to relate TD errors to position errors; l.e. , 
the angle of the LOP passing through the receiver position Is broken 
down to cartesian components and the TD eiror is used as a weight Li\g 
factor along with the LOP components to calculate the new position. 
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Figure 3-3. By Movif>g the Receiver Posirion (R) Closer to the Mosfer (K*) 
Along the Moster-Receiver Baseline, the Assumed TD May Be 
Mode to More Nearly Correspond to the Received TD 
Provided the Initial TO was too high. 


- 26 - 



Fi{jur« 3~4. A Mefhod of Convorting TD» to Position b/ Minimizing the 
Gxnporison Error by Moving the Assumed Position Closer 
to or Farther from the Moster. 
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This procuss is iterated uatil the TDs bused on the adjusted position 
agree with the measured TOs within an acceptable range. 

The mathematics of this algoclthm will be more fully developed 
at thl point. Keferrlng to Figure y-b, tlie tangent to the hy))erbuUo 
LOF p<»saing through the receiver position is: 


f 


SM 


2 


(3-1) 


where and are the bearings , from true north > of tlie geodesic 
arcs from the receiver to the secondary and master stations, respec- 
tively. The tangent to the hyperbola also bisects the angle between 
the two arcs to the Loran-C stations, as shown in Figure 3-5. A rela- 
tionship between TD error and position error is found by considering a 
vector whose magnitude is: 


1 



(3-2) 


The north and east components of | V | are [22] : 

a ■ ^ (costg - cosS'j^) - -Vsinfgi^ 


Y - ^ (sintg - slnf^) - Vcosfgj^ 


(3-3) 


The TO errors are calculated as follows: 


6T 


A 


‘TA 


OA 


“ ^TB ■ ^OB 


(3-4) 


where T^^ and T^^^ are the TDs corresponding to the assumed position, 
and Tq. and T^^ are the measured TDs (a position fix requires the use 
of at least two secondary stations, resulting In two time difference 
measurements). TTie new estimate of position is: 
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Figure 3-5, G#om0rry Mcuter (M), SoconJory (A) ond Keccivcr (K) 
Polifion. The fooyuiU to the LOP bikoct^ the umjle 
formed between the geodesic ore from the receiver tu the 
mo&ter and secondory. 
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‘2 ■ *1 - ‘Vi> 


“2 • “1 ^ - ”a - "«“2> 


wheru 


D ■ ^*^^2 "* ^2^1 * ^ '^SMA * ( 3*u) 

and Che eubecrlpte I and 2 refer Co quant It ies associated with secon- 
dary atationa A and B» respectively. 

Analysis of this algorithm showed that it involves tairly 
extensive calculations; as many iterations are required* the Tl) pre- 
diction algorithm is called many times and can cause a severe time and 
memory penalty. 

For use in a navigation unit* once a |X>sition fix t\as been 
determined I succeeding fixes will he close to the preceding ones* a 
fact which can simplify posit Ian fix calculat Lons , Htuxy iterative 
algorithms cannot ^lake use of this fact* since their first correction 
may overs!\oot the actual position at\d nkiy be along a direction which 
is not coincident with the ivivlgation track. 

A computer program developed by the Haval t>ceanographlc 
Office 1231 uses a technique whereby a set of Loran-C Tl)s may bt* con- 
verted to latitude and longitude if a starting |X)int close to the 
actual position is known. This process uses the baseline lengtlis t*> 
the stations (expressed in microseconds) to set iip a homogonoi>us si'i 
of transcendental equations which are then solved tor unknown correc- 
tions applied to the initial position. If the initial }X>sition Is 
more than 10 to 20 kilometers from Che actual jx>sition, cyclic ambi- 
guities arise in the solution to the set of equations, thus a series 
of iterations on the initial position may be necessary before the 
actual position may be determined. This program will converge 
directly on the correct position, thus eliminating the problem of 
having a large radius of converge* common with Che previously 
discussed cases. Although this program Ivas proven to he very accurate 
for use on an IBM 370 124), a considerable number of calculations are 
necessary to obtain the solution. For use It a small computer, this 
would result in a large amount of memory for the cmlculat ions. The 
solution is faster than the previous case; however* since only one 
iteration is necessary once an initial point close to the actual 
receiver position has been determined. 

Several other iterative-type solutions are |>ossible and leave 
been investigated. One set of conversion methods in [xtrtlcular has 
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been developed which solves the TD^to-poslt ion on the transverse 
Mercator xy grid and are fairly simple [25]. However, to use thin, a 
large set of conversion equations must be solved to convert from the 
UTM (Universal Transverse Mercator) coordinates to the standard 
lat/long coordinates^ which Involves the solution of many transcenden- 
tal equations. 

A non-iterative, explicit solution to the TD-to-posit ion 
problem has often been envisioned as an ideal substitute to the above 
algorithms. Explicit solutions are complicated by having to deal with 
the non-spherlcal nature of the earth and also witli the non-constant 
propagation properties of the Loran-C signals, as discussed in Chapter 
II. These effects can be Introduced as local correction quantities in 
the iterative solution once the iterated position becomes close to the 
actual position. However, with explicit solutions, these corrections 
must take the form of a general model v;hich represents a large area. 
Thus, an explicit solution can be simple if precise accuracy is not 
necessary. 

Non-iCerat ive solutions may take several different forms, 
however, most form a relation between the received time differences 
and the corresponding distances from the receiver to the Loran-C 
stations. From here, these distances and angles between baseline arcs 
are used to form a final expression from which the receiver position 
is found [26]. 

An excellent method for converting TDs to |X)8ltloa coordinates 
In the form of latitude and longitude was presented by Razin [27] and 
Is presented here. Referring to Figure 3-6, the spherical arc angles 
to the secondary stations are: 
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where : 

- -R - ■fcx> - 0«x 

■■y - T - ■'ey) - % 

(3-8) 

and V « v^/n, the velocity of light divided by the index of refraction 
over the surface, R is the radius of the sphere in Figure 3-6, and 
T are the measured time differences received for stations X and Y, 
respectively, and T^^ and T^^ are the secondary coding delays for sta- 
tions X and Y, respectively. Notice that in this case, the measured 
time differences are converted into spherical arc angles for the 
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Figure 3-*6. Geometry On a Sphere of a Loron-C Complex Con^iktin<j 
of AAoster (M), SeCortdory-X (X), Secondory-Y (Y) ond 
the Receiver (S). The 6'$ ore arc ongiei meosured at the 
center of the sphere ond the other angles ore sphcricol 



r«c«lver buitellnes to the Loran-C statlone. Now taking the coKlne of 
equation (3-7): 


- sinO^gSlnP^ 


co»6y„ " <^o»0^coBPy - »inO||jg«lnPy 


(3-9) 


Proa the properties of arc lengths on a sphere, these equations nay 
also be expressed as: 


cosOjj^ cosO^cosO,^^ + slnOjj^sinO|^cos(J^ 


cosOyj^ "" slnO„ySlnQ„,^cosfiy 


Equating equations (3-9) and (3-10) results In 


tanO 


cosP - cos0_„ 

X IkX 


«8 slnP^ + slnO_„co80, 
X nx X 


Fr<xi Figure 3-6 . 


(3-10) 


(3-11) 


cosSy ■ co 8(K-0^) • cosKcosS^+slaKsinS^ 


(3-12) 


and 3^ aay be solved in terns of the known quant L t ten : 


cos6 - 
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u^ij^ + rUj U2 “ 
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^2 ^ 1 '' ^^2 
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“l -^2 


(3-13) 


tfhere : 


Uj - a^^cosK-ay 
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X mx 
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y 


cosP “ cosO 

y ^y 

s tnO 


my 


b 

X 


siuP 


sinO 


X 

mx 


(3-U) 


Tht» angle 0 . is found then by: 

ms 


0 

ms 


tau ^ T — 


0080^ 


(3-l->) 


ind equation (3-7) may be used to find 0 . and 0 , 

xs y s 


To find tht? latitude and longitude of the prosenl iHjsltton, ttu- 
spherical arc angles Just calculated must he converted to spherical 
latitude and longitude. This process Is rather compl Lcated , but the 
results are a scries of constants multiplied by the arc angles. Thi’sv* 
constants need bo computed once for each set of slat ions used u\d nviy 
then be stored In the navigation computer. Let and he tle^ 

latitudes of the master, secondary-x and secondar*y-y and 0 , ^0 and 
0 be the longitudes of the master, sccondary-x and secondary-y. T\w 
receiver’s position Is denoted by 4> anil 0. The following relations o. 
a sphere wA*re developed by fell 128) : 


cos^ sinO (cos4>sln0) -f cos| cosQ (cos4»cos0) + slni^ sln<^ - cosO 

( i-U>) 

co8^y!,laOy(co8^slnO) + co.^yCosOyCcus^icoHO) + ttlii4,ySla<(. - «-usOy^ 

(3-17) 


Oyji and from equa*:lon8 

0-15). 

+ C3cos0y^ 

0-19) 

+ VosOyj, 

(3-20) 

+ C^cosOyg 

(3-21) 
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cotf4p^silii0 (cosAslnO) + cos4 cosO (cosiicosO) + siuO siaii. - cosu 
® m m ^ m ms 

OIM) 

Theiit nay btf treated a« three llaear equations In tiireo unknowns by 
letting g-co«4»sin0, f^cos^cosO and h*sin4i, line solutions of < liesi* 
equations are found by using Ox 
(3-16) and (3-17). These are: 

f - ^ C^cos 

g - C^cos 

The constants are: 

C, - (cos6 slnO_sln(>„ - sln^ co««l>„sinO,,)/D 
1 m y n 'y y 

C 2 • (sin^j^coH^ySinOy - co84>j^sin0j^8ia4i^,)/ 1) 

C« - (co84„slnO„Kln4i„ - cos* sliiO„)/D 

J XXX X m III 

- (sin^^coBi^yCosOy - co8v(#^cos0^sin4»y)/L) 

Cc ■ (co8^ cosO sin^ - sin(^ cosO cosO )/D 
j X X. y X y y 

Cr " (slnA cosA cosO - cosd> cosO slnd> )/Ly 
u x ni ra X x’^m 

- (cos A cosA^)(co80^sin0,, - sinO cosO„) / 1) 

/ to y my roy 

Cg - (co84i^co8^y)(sln0j^cos0y - cosOj^s UiOy)/U 

C_ ■ (cosA cosA )(cosO siiiO - slnO cokO )/U 
7 x m xm xra 


0 - 22 ) 


D - ‘*irt|>^cof.^^co8(>y{coB0|j|Sln0y - slnO^cosOy) 

- Bln* co8^ co8(> (cosO slnO - slnO co.sO ) 
ID X y X y X y 

+ sin* cos* cos* (cos 0 sinO - sinO cok(j ) 
y X n X m x n 


(3-2'J) 
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Then Che latitude and longitude of present posit ion is; 

4» » sin (3-24) 

0 - tan"^ (3-2‘j) 


The conversion of the received TDs to position us ini’ this 
iLetnod is illustrated in Figure 3-7. This algorithm uses a spherical 
model of the earth to ttnd the lat/long of present position. By 
extending the treatment of converting the spherical arc angles to tl\c 
corresponding geocentric coordinates, an elliptical mode L of tlie earth 
may be incorporated. 

Using the Clarke 1666 spheroid, the major radius of tlie earth 
is R and the minor radius is r. Let with latitude and longltu- 

tude OQt be approximately in the center ot Che coverage area of the 
master and two secondary stations. Hien the center of t'ne sphere,* 
osculating at Pq (see Figure 3-8) is: 


- (R - i/R^sin^4>^ + r^cos^ij)^) cos<^^^slnO^ 


« (R - -F r 


2 0 ^ ~ 

-F r cos*'4>q) cosO(jCos4>^^ 


- (r - -5 + r^cos^ 




❖o • tairU-p tan^)^) 


( 3 - 26 ) 

( 3 - 27 ) 


(Geodetic, or map latitude, is denoted by (>q while geocentric latitude 

by Also let p be a point on the spheroid with latitude 4) and 

longicude 0. 

Its Cartesian coordinates are: 

X - Rcos^binO 

y - Rcos^cosO (3-28) 
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Figure 3’'/. OepicHon or Explicit TO-to-Positior> Algorithm. The dotted 
lines represent corrections applied Tor propogotion effects 
'^'><3 ellipticul earth corrections. From discussion by Rozin. 
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Figure 3-*8. Mapping of Point P On the Clarke Lliipsoid 
to An Osculating Sphere Whose Surfoce is 
Tangent to the Point P. 


iiusaiffieifai 




OF POOR QUAUTY 


Let p’ be the liiMgo of p under the nkii^ trow the sphorviUi to tho sphorv** 
Then the coordlnetea of p' ere: 


*x) ^ r^. (y-') ^ r^,(E-Z) 

1 . I, 

(3-30) 

((x-X)>‘ + (y-f)^ + 

( )-3l) 

• r\o«' 

(3-K-) 


r 


T\\e latitude and li'ngltude of p* are: 



0 * 


ta.r‘[ 


x-X >1 

y JY J 


( n) 


To lmplenH?iit thla, the above procedure Ih doiK' tor each ol Uic 
a tat Lone. In other ^rdH, letting 0 ^* 0 ,^ and Oy be the 

latitude and longitude of the laaater and two secondarien, then the 
image of these on the sphere are calculated giving 

0 ^ and tV * These are then used to calculate the constants througti 
as nefore. 


With these nine constants and the calcuUUed arc angles, the 
latitude and longitude must be nsipped back onto the spheroid, ’Pu* 
result is; 

-I .K + 

- taa lC|^,ln 0 g, g'+ ( 3 -JA) 
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( i~y>) 


The other coobtaatMg through nre ; 



^13 ■ 

^ . R (3-30) 

7 

Note that the computation ot the C-constHuts Is done *once for 
each triad of Loran-C stations and In the navlgatl'^i computer; only 
equations (3-7) through (3^20) are computed, and the C-constants are 
then used as mult Ipl lent ive ^>arameters« A FX)RTKAN-LV program lmplO“ 
mentlng the non-1 terat Ive TU-to-posit Ion procedure Is shown In 
Appendix II. A FOHTRAN-IV program which calculates the C-constants 
using equations (3-21) through (3-33) Is given in Appendix IV. U can 
be seen that the TD-to-poslt Ion program is fairly short and does not 
Involve very many complex calculations. Figure 3-9 shows a breakdown 
of the mathematical operations used, 

D. Distance^ ajuI Corapujt^at lons . Having presenL 

position described In terms of latitude and longitude Is of limited 
value for the pilot since the pilot will often ru'cd to refer to a 
chart or map to locate the present i>osltlon relative to geographic 
features or other navigation aids. It Is more helpful to have the 
present |X)8ltlon given in terms of the range, or distance, and the 
ailmuth, or bearing angle, to or from a waypoint, which could be the 
departure airport, the destination airport or some other reference 
point along the aircraft’s route. The pilot may wish to use sever,il 
different waypoints and these would in turn be stored in the ivivlga- 
tlon computer. Then, as the flight progresses, the computer would use 
the currently selected waypoint and the iH?west estimate of position to 
compute the range and azimuth. A way of computing this was presented 
In Chapter 11 for the development of computing Loran-C TDs given a 
known position. An alternate method Is presented liere which Is 
simpler and better suited for these range and azimuth (also called 
rho/theta) calculations. Referring to Figure 3-10, the coordinates of 
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Figure 3*9. Breakdown of Morhemotical Operotioru and Voriobies 
Lead in Explicit Coordinate Converter Progrofn. 
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[ Figure 3«10. Geometry on Sphere Lhed for Range ond ikioring Angl' 

I to Waypoint from Present Position Calculations. 
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prutiuaC |)OttlClon are denoted by and and the coordinates of the 
waypoint are denoted by and X^ for the latitude and longitude, 
reapectlvely. The bearing angle to the waypoint is [29] : 


♦ru " 


-I 


♦ r ♦u 

(xr - x„) co8(I5_-::«) 
♦w “ ♦r 


(3-37) 


aad the range to the waypoint is: 


Srd 


- R (Xr - X„)cos(l«— + (♦w - 


1/2 


(3-38) 


where R is the radius of curvature of the earth. This rnay be coioputcMl 
for the oU.d-region of the coverage for the particular Loran-C stations 
used and stored as a constant. After choosing the midpoint of the 
coverage region, R is found by: 


2 2 2 2 
a ♦mn cos 

R 

b (3-39) 

where a and b are the semi-major and semi-minor radii of the reference 
ellipsoid and is the latitude of the mid-point. 

E. Pro gram f or C oordinate Conve rsions . Table 3-1 lists trade- 
off considerat ions for the various TD-to-position methods. The non- 
iterative TD-to-position procedure outlined In Section C was found to 
be the most suitable of any for application to a microprocessor-based 
coordinate converter. As already mentioned, the most complex calcula- 
tions, Involving the C-constants, may be carried out on a central com- 
puter. Then the resulting C-constants are transferred to the 
microprocessor system. For any given triad of Loran-C stations, the 
same C-constants are used regardless of the aircraft's position. 

The ellipsoidal model gives better results than the spherical 
model, as might be expected. The degree to which the ellipsoidal model 
does better depends on the choice of Pq in equation (3-25). This 
should be chosen to be in the center of the coverage region of the 
triad. This Is most simply dene by averaging the Latitudes and longi- 
tudes of the master and two secondary stations. This average latitude 
and longitude is then used for Pq* 

The direct TD-to-position algorithm described in Section C was 
tested for accuracy and to determine if it exhibits any singularities, 
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TU-to-l‘obit lort 
Method 


Advtf nlH|*ee 


UlbadvaiU .i^;i*s 


Iterative Method 
of adjuetiac i>oeition 
cloeor to or Carihor 
f ru«a tAiMlyr. 


Simplo to prograio 


Could rc<|uirc 
iar^c oumhoi* ol 
1 tor.it ions lo 
Convur4;i: on result • 
TImo conb»nnln^;. 
Lar^u radiiis ol 
coi^ ve r^;o«iCv* • 


Iterative nethod of 
edjuetlog i>oeltlon 
beeed on else of 
TO coMperieon error. 


Requires fewer iter* 
etione than above. 
SmH radius of con- 
vergence. Possible 
to insert local 
propagation correc- 
tions for better 
accuracy . 


Tine consuming - 
TU prediction 
rout Lne c.i I le«l 
repeatedly. Re- 
quires More uieiaory 
than above. 

Possible Lo have 
|H)inis ol singu- 
larity in solution. 


Direct one* its rat ion 
solution for assuMod 
position close to 
actual position. 


Provides high accu- 
racy. Paster than 
above two. Result 
is given with pre- 
dictable oaount and 
direct ion of i>os- 
slble error froa 
true !>osition. 


Much larger •is*ia«»ry 
requi re4k.*ni s than 
above two. Ke(|uires 
iterations or dead- 
reckoning it assumed 
position is i.ir I rum 
actual iHjsUion. 


Direct closed forui 
solut ion. 


TD- to-pos I t ion cal- 
culations are laucli 
sluiplor than above 
three; done once. 
Highly accurate. 
Predictable error, 
if present. Modest 
Beiiory storage 
requi reMont s . 

Faster than above 
three . 


Spheric.il c.nglc tv» 
geocentric |K)sition 
caicul.ilion coi^st.inis 
need t*) bo computod 
of I line. More dUl i- 
cult than with 1 ifsi 
two lo Insert propa- 
gation corrections. 


Table 3-i. Tradsoft Cons ids rat ions lor Various 
TD-to-Pos 1 L ion Methods. 














or puialM ^i:re convergence wae (>oor. The results of UiLh LeiiL .in; 
ehown in Figure 3*11 » which ehowa randonly aclecCed polntM along with 
their CO r re g ponding TDs, which were calculated using a program devel- 
oped by the Naval Oceanographic Ofllce» and the calculated position 
determined by the direct algorithm using the program in A|)p€ndlK lU 
TVo tests were conducted; one using a spherical approxiiMt ion of the 
earth and the other using an ellipsoidal approximation of the earth* 
As can be seen, the ellipsoidal version gives more accurate solutions 
and the solutions tend to show smooth deviations from the correct 
values^ as shown In Figure 3-12* Also, as evident by standard 
deviations, the spherical approximation shows sharper deviations from 
the normal errors expected* 

It is quite evident that the ellipsoidal loodi*! gives better 
performance than the spherical model* The two models were run using 
FORTKAN-IV programs of the type shorn in Appendix II on an IBM 
370/138* The programs used approximately 1300 byvtH of memory and 
executed in 0*02 seconds* 
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Figure 3”11. Test Results of Lxpllclc Coordinate Converter Fr'>:{rau 
Using Both Spherical and Ellipsoidal Earth Hodeis. 
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Figure 3-11. Contlaued. 


- 47 - 



ORIGINAL PAGE 13 
OF POOR QUALITY 


C«nl«r of 
Covoroge 


G«n«ral RoloHonihip of Convoriion trron to Di»tonc« from Central 
Point. The top number in eoch circle is the sphericol mockfl error, 
and the bottom number is the ellipsoidol model error. 
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This chapter will detail the process of running the direct TIK 
co*poaltlon program described in the last chapter on a micrcprocessor- 
based system, k description of the microprocessor system to be used 
will be given along with the necessary additional hardware required. 
Also, special software problems unique to this process will be 
discussod. Pinallyi results will be given showing how wll the coor- 
dinate conversion program works and Yu)w well it can he used with a 
Loran-C software-based sensor processor receiver in a general aviation 
environment. 

Host popular micJ^ processors on the market contain an accumul.i- 
tion register (accumulator) which is eight bits wid«*. Ail data pro- 
cessing, then, takes place on eight-bit words. These prcn'cssors 
contain* other registers; for instance, most have a l6-bi register 
which drives the processor's address bus to fetch instruct ions and 
data, there are us*ially several registers internal to the processor 
which can be uaed for scratch pad purposes in calculations or as 
indices to an address, and there is usually a register Cor Lnoicatlng 
the status of the processor at any given time. Most processors are 
not suited to doing arithmetic work such as miilti-prec Islon addition, 
subtraction, multiplication and division when only eight-bit words can 
be handled at one time. Because of the range of numbers involved in 
the calculations in the coordinate converter, It is necessary to use 
f loat ing-point numbers; arithmetic computations on these can be rather 
complex. 


A microprocessor which has been used previously in l^oran-C 
development work al Ohio University is the HOS Technology H0S6302; the 
Internal architecture of this microprocessor is shown in Figure 4-1. 
This microprocessor con«iate cf one accumulator where arithmetic and 
logical reaults are collected, two eight-bit index regfsters, one 
eight-bit stack oointer, a six-bit flag or status register and a 
program counter [30]. Figure 4-2 is a list of instructions n\d 
addressing modes for the 63^2 [31J. The softwaie development support 
for this system includes a symbolic assembler [32] and a simulator 
botli of which refolds on the Ib:i 370/13tt at Ohio University. Host of 
the microprocessor work done ar Ohio University %#as done on a Super 
Jolt microcomputer, shown schematic nil y in Figure 4-3. 

l>€slgn C ons ide rat ions. The coordinate converter algorithm 
described in the iMt cYTa^pter ms written in POUTKAN-IV for testing 
and evaluation. In order to run the coordinate converter on the Super 
Jolt, it is necaaaary to convert the FORTRAN' IV program statements to 
assembly language atatementa which the 6302 can then run after the 
statements have been assembled. Tlie main problem to he solved la 
doing this is the mathematical operations which must be perforiaed. 

A table is set aside In memory wtiich ccntalne various numbers used as 
constants end also numbers uaed for lntermedl«‘ite calculations. As was 
menrioned, it is necessary to use a floating-point format for repre- 


























'oressina 



















ttcuting the immh^irs since the numerical range Is .(p|>roxlmate I y 10 
10 • approxl«HCely Il5 orders of loagaltude. The software to do 
floating-point arithmetic on a microprocessor uses a large .imouiu of 
memory and can take several milliseconds to do a f lost Ing-po Ini 
mul t Ipl Icat Ion If the dlflerence hietwcen tlie r.inge ol numbers is 
large. 

Tlie /unoiint of memory in most small mlcroprocesst>r systems Is 
limited^ especially for applications requiring much aritlimetio pro- 
cessing and interfacing to other <Ievlces, Tlie microprocessor system 
used at Ohio University for Loran-C develoi>«ent work Imk flv#* kilo- 
bytes of read/write memory available and provisions for >idding more 
memory, either read/write or oermaaent type. To cut down si/.e and 
cost, It Is desirable to place the coordinate converter control 
program in a 2048 byte erasable, read'^only memory (EPKOM) and use 
approximately 500 bytes, or less, of read/wrlte memory for scratch pad 
calculations, storage of computed results, and transfers of data to 
and from the scnsor-prucessor. This represents a practical require- 
ment from the sense of designing and building such a system. 

As piirt of the project, a peripheral mathematics processor was 
used to do the necessary floating-point calculations. This device is 
the Ain95ll by Advanced Micro Devices and Is shown In Figure 4-4. Tills 
device Is designed to be used on nilcroprocessor systems which use an 
eight-bit data bus. It can handle l6-blt and 32-blt fixed |>oLnL 
arithmetic and 32-blt floating point arithmetic. Figure 4-5 Is a list 
of instructions for the Afn95liA. 

The input to the coordinate converter comes from a 1/iran-C sen- 
sor routine which runs in the same microprocessor system 133). Tlie 
output of the coordinate converter consists of the computed lat/long 
for the received TDs and the range and bearing to the selected 
waypoint. These values are dlspl*iy^d on a video display unit con- 
nected to the microprocessor systeii ^4]. These computed valmo;, 
aloii^ with the ntceived TDs are also recorded on a data n!Cor»liii); unit 
which is interfaced to the microprocessor system. The overall scheiiu* 
of the microcomputer system used for the coordinate converter is shown 
in Figure 4-6. 

B. Syste m fo r Coordina te Coni^rsion. In order to Interface 
the Aiii95llA to’ the JoTt microcomputer system, it was neccHsary to add 
additional hardware to tl-c Jolt to handle device selection and data 
transfer, and also to write a set of subroutines in assembly language 
to liandle the transfer of data in user memory to the 9511 and vice 
versa, to conunand the 9511 to perform a function and to check the sta- 
tus of the 9511. Since the 9511 was designed to Interface mainly with 
the 8080, 8085 and /.80 family of microprocessors, wlilch use a dif 
ferent method of handling data transfers and monitoring of [lerlpher.il 
devices, an H6820 peripheral interface adapter (PIA) was used for the 
hand-shaking. The method of doing this is illustrated in Figure 4-7. 
The M6820 consists of two eight-bit [xirts and several other registers 
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Figure 4-7. Schematic Diogrom of Microcomputer Intcrfoce to VSI I 
Oeiigr^ed for TO-to-po^ition Computation S/&tcm. 
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iisn! lo to |vri)>hcr.il ilevlcfs .is shown in Kl;;nr»* l 

Slili* \ of the PIA irt nsocl to transfer eight hits ol •lata tin* 

9Ml and the micriiprocossor system. Various |>orts on side li are nse»l 
to set and read I/O lines on the 9^11 for rending dita, writing d.ita, 
selecting Internal registers In tlu» ‘ibll, etc. 'Hu* interrupt regls- 
t%*rs on the Hf>H20 are usiul by tlie 9*> I I to Indicate when llie 9Sll Ivik 
made valid data available to the microprocessor and when the 9bll li is 
completed eaecntlon of a command. 

Tlu‘ next part of the Interf. icing to the 9511 Is tlie de ve lt>pmenl 
of a standard set of progr.uiis ( suhrout Ines) which an^ used to Ini- 
tialise the lurdware Interface and the 9511, to write, or load a 
slng.lc f loa t Ing-po Int (32-hlt) number, to read or uolt>ad *i float log, - 
oolnt (32-blt) number, to send an elglil-blt word to the 9511 repri*- 
sentlng a command to bo executed, and to road the 9511 elght-blt 
status register. Figure 4-9 gives a summary of the necessary proce- 
dures to do the Items Hated above. It Is desired to make these 
subroutines complete to the ixilnt where the iiwiln program can specify a 
pointer to a ivirtlcular number or specify u )>artLcular operation to Ik* 
executed and then call the appropriate subroutine to do the function 
and return when the operation is completed. To this end, an area of 
storage In read/write memory Is reserved for storing all numbers which 
will eventually be i^assed to and from the 9511. The biise address 
(beginning of table address) is stored at an address to wtileh the 
subroutines have common access, then all tlie numbers are referenced by 
a single elght-blt offset from the base address. Since e.ich number 
occupies 32 bits, or four bytes. In memory, 64 numbers can be stored 
by tills procedure. If more than 64 numbers must be stored, tlien a 
sec«>nd nuinbi^r table and a different base address must lie sol up. 

Figure 4-10 shows a set of flow diagrams to Illustrate the 
logic flow of a set of standard subroutines to do Che processes pre- 
viously discussed. Tlie first, called **PINI\*" initialises the M6M2n 
Pl/V, which consists of setting the data |X>rts to either Inputs or out- 
puts, as required, aiiil setting up the Interrupt registers, then the 
control lines to the 9511 are deselected so chat the 9511 will not 
respond to Invalid data during the In i t ial I xat Lon iK^rlod. Also 
Included In this subroutine Is a section which moves a set constants 
stored In the |>er»uanent memory to the read/write numher table anil Inl- 
tlall/.es in ari*a of ixonory used for vldet) display work. 

Subroutine **PUSIl‘* is used L*» copy a four-byte number from 
reid/wrlCe memory onto the stack of the 9511. The Y-reglster contains 
the eight-bit offset from the base address of the number table for the 
number to be copied. In this case, four Individual transfers are 
made, since the number is 32-blts long and only eight bits are trans- 
ferred at a time. The nuiaber which is copied appears at the top of 
the stack In the 9511 and any other numbers stored on the stack arc- 
moved further down in the stack by one location. llie number stored at 
the bottom of stack Is lost. 




- 58 - 




















ORIGINAL PAGE IS 
OF POOR QUALITY 

Upor.iiul ciUry 

1, I'l.it'o Ic.iKl b 1 iC4iiil UyLi* on il.il.i Un.. 

2. Scl i:/d low. 

i. ScL low. 

4. Sec WK low. PiAUbe will tH:coiue low. 

*). Piiub^ relarrib lilgS. 

0. Set WK hl^h-data, C/I)\ CS inay Chen 

/. Repeal above bix Lo write entire 32" bit number. 

Dat.k removal 

1. Sot C/ i) low. 

2. Sot CS low. 

3. Set KD low. Paube becomes low. 

4. Pause returns hljjh - datii is then available^ must sli;nil i- 
cant b yte lirs^._ 

5. Set KD hlijb - CS, K/D, data iiwiy tlion clkanj^e. 

6. Repeat above live to read entire 32-blt nuinl>er. 

Conuiund initiation 

I . Place ^oiomand byte on data bus. 

2. Set C/D lili»h. 

3. bxccute steps 3-b ot operand entry. 

4. ^ND output remains low until coiaman * completes. 

Status read _ 

1. Set C/D 

2. Kxecute steps 2"S ol dat.i removal. 

3. Status byte is available on bus. 


Figure 4-9. Summary of Sfep& to Communicate to 9511. 
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Figure 4-10. Cont’oued. 












Suhri>uCliic **POP** l8 ubCiJ Lo do the opposLle ol the •ihove ; *i 
32~blt numbor in removed from the top of btack in the 9^11 .ind 
In btorage. The addrebb where the nomher lb stored Lh |iolnted l<» by 
the Y-regl^ter offset relative to the habc addrebs. After tills opera- 
tion lb coieplote, the nuniher that was at the top of stack is lost nul 
any other namhers in the stack are incremented up by one |K>sition» 

This subroutine incorporates a test of the [lause output of the 9bll 
before loading datK over the data bus. In Home instances, It may take 
the 9511 several microseconds to present the data; when the data Is 
available, an interrupt bit in tlie PIA is set. This Is not necess.iry 
Cor writing data, however. 

Subroutine **CMND'* is used to couiuiand the 9511 to ix^rCorm .i 
given function, the function code is contained in the nccniuul.i tor , 

After the 9511 has ^iccepted the command, it immediately b«fglns esecn- 
tlon of that command. During this time, its status register will 
indicate that it is busy; when the command l^s finished execution, the 
status line will set an interrupt bit in the PIA. The Hubroutine 
loops until this hit is set before returning to the imln progr^n; thus, 
on return, the function will be completed and the isaln program need 
not check for complett'a. 

As |sirt of the **CMNU** subroutine, a check Is made of the siiLus 
register to determine the final outcome of the completed command. 

This may consist of an error In processing h4*cause of illvislon by 
zero, an input number out of range, etc. The return code is reteruod 
in the accumulator for further processing hy the miln program. 

Other subroutines are necessary Cor the conversion of blurry 
coded declisal (BCD) input data representing the time dlfteren.i's Lu 
the internal floating-point format used by the 9511 and frots ’*• '-ting- 
point to BCD for displaying the final results. One of the c^.km-aJs 
available with the 9511 is one for converting a fixed-point (integer) 
number to floating-point. The fixed-point number mos;. I>e in binary 
and it is not possible to represent a fraction using this command. 

The time differences from the receiver are in base 10 in units of 
microseconds with a single digit representing a fraction of a microse- 
cond, for example, A2594.2. The two lime differences are stored at a 
certain location in memory by the receiver and are In packed format; 
i.e., two digits occupy one byte of mi?mory. Before the TD-to-pos it ion 
calculations may be performed, these two TUs must be converted to 
binary, then to floating point. In tlie microcomputer program, this is 
done by treating each TO without regard to the the decimal point; the 
actual binary TD is multiplied by 10, the correct Tl) is obtained by 
multiplying by 10*' to get them in microseconds. 

Pie piocess of converting a base 10 number Lo binary is 
llluHtrated In Klgure 4-11. This is done by starling with tlie most 
significant digit (MSD) of the base 10 number and adding it to a 
register then |>erforinlng a binary mul t Ipl Lcat ion by 10. in this way, 
each digit position of the base 10 number is converted to binary by 
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Lite mult Lpl Lcat Lon and added to a pari Lai sum to form llie final 
answer • 


The binary raul t ipllcat ion Is performed by a process called 
liooth's algorlthn [36] • This is also applicable for binary division, 
as used at the end of the lalcrocomputer coordinate converter program. 

In this algorithm, the multiplicand is stored in successive locations 
in memory and the multiplier is stored in the accumulator. The 
multiplier (in binary) represents a string of zeros and ones; the 
number of iterations through the multiplication loop Is set by the 
length of the accumulator, in this case, eight. With each Iteration, 
a residue table (Initially set to zero) Is rotated left, the bit 
shifted out of the left, reenters on the right. Also, the hits In Llie 
accumulator are shifted left, %dien a one Is shifted out, the multipli- 
cand is added to the residue, if a zero is shifted out, no adding 
occurs. This is illustrated in Figure 4-12. This process is very 
similar to the way multiplication would be li*indled with decimal 
numbers; as the residue is shifted left, It is effectively Increased 
by a power of two. When a one is encountered, the multiplicand Is 
then added to the partial sum. 

When the coordinate conversion process is completed, the 
resulting latitude and longitude of present position and tlie range and 
bearing to the waypoint are in the binary floating-point format and 
must tiien be converted to BCD for display and use by the l/>ran-C sen- 
sor program if they are to be stored on tape oi* processed furtlier. For 
the latitude 4ind longitude, the coordinates must be converted first 
from radians to degrees, then broken down to the standard degrees, 
minutes, seconds format* This is done by taking the Integer (virt of 
tlie coordinate ami storing It as the degree i^art, then subtracting tin* 
integer from Che original coordinate and multiplying by 60 to obtain 
the raLUutes part. Doing this a second time yields the seconds |>arL. 
Most of this work is done by using Ciie floaC-to Fix and flx-Lo-float 
commands In tlie 9511. 

The binary to BCD conversion is done by taking tlie binary 
number and |>erformlng a binary division by ten. This yields a 
quotient and a remainder, the remainder Is the cor respond Lng decimal 
digit of the new base 10 number, working from LSI) to MSI). The 
quotient Is then divided by 10 again to yield the next decimal digit. 
This Is continued for a number of iterations corresponding to the 
number of bits in the original base 2 number. The subroutines which 
perform these functions appear near the end of the program listing in 
Appendix IV. 

The Loran-C coordinate converter program was designed to run on 
the Super Jolt microcomputer a ng with another control program which 
operated softwaie aided loops that lock onto the Lorau-C pulses and 
compute the correKponding time differences. The executive program in 
this also operates Interfacing to a video display and a data recorder. 
Certain addresses in page zero (0000 hex to OOFF hex) of the MCS6502 
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nieinory bcUcmc aro reserved lor d.ita coinmun Lea t ion between c*n -.ni-C 
sensor program and the coordinate converter proy^rain. 11 eMe ...» ».ndi- 
Gated on the first few pages of the program listing in Appendix IV. 
These locations are used to pass the two time d i f fereiui*s to the c«>«)r- 
dlnate converter and to receive the subse<jiient latitude, longltiuU*, 
ranj'i* .ind bearing. All of tliese are stored «is packed IICI) nuiuhiTs. 

The complete mlc roprocessor progr.nn for tiie coordinate conv«.M‘- 
sion Is shown In Appendix V and occupies 2043 bytes of KPKOM. Hu* 
progr.iin also uses one^half kilobytes of read/write wmory for scimIHi 
calciila t ions. Figure 4-13 shows a photograph of the control hardw.uc 
including the Am93ll, the 6820 PIA, a sp.in? KPKOM and associate*! 
inter facing components. This board plugs int*) tiie I^r/in-C developmen- 
tal r»-.celver, ihown in Figure 4-14. 'PU s unit was LesLe* noth on Liu* 
bench and in several flight checks, which are reported in the iu*xL 
chapter. 


As an aid to preparing the coordinate conver*:er software, .i 
service progr.iin was written for use on the I KM 370 to allow b.ise in 
numbers to be ^yped in from the user’s terminal, and Llitui converted L*> 
the equivalent binary floating-point mimber for the 9311 and «lispl.iy*.»tl 
on the us*.»r’s terminal. The program listing for this ippe.irs in 
Appendix VI. 
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V. TESTS ON MICROCOMPUTER SYSTEM OF POOR QUALITY 

Details of the developaent of the Loruii-C cuordiaate coiwerler 
and the microprocessor system have been given. Foliowing is a verifi- 
cation of the performance capabilities of the coordinate converter 
system. This consists of tests run to determine the accuracy 
the process and the suitability of the system to used in an aircraft 
environment for navigation. 

The microcomputer coordinate converter was placed on one 27 in 
erasable » programmable, read-only memory (EPROM) chip. This EPROM luis 
capacity for 2048 bytes of data; the coordinate converter occupies 
2045 bytes. With the Am9SllA running at a clock speed of I MHz, a 
complete conversion requires 0.18 second to complete. (Ttie maximum 
clock speed on the 9511 is 2 MHz and data transfer between the 9511 
and the central processor is asynchronous, thus computation throughput 
could be doubled by running the 9511 at 2 MHz while keeping the 6502 
central processor at 1 MHz.) The total memory requirements for the 
coordinate converter are: 2045 bytes of permanent (EPROM) storage for 
the control program and constant numbers, 37 bytes of page-zero 
storage for temporary variables and flags and communication to the 
Loran-C sensor program, and 120 bytes of general read/write storage 
for results of intermediate calculations. 

Testing of the coordinate converter consisted of checking 
repeatability, ability to compute correct position given the time 
differences of a known point, and a check for any singularities or ten- 
dencies to converge on the wrong result. Checking of repeatability 
consisted of using the same time difference readings as input and 
determining if the coordinate converter gave the same answer each 
time. In all cases, given a particular TD reading, the coordiivite 
converter yielded the same answer regardless of when the test v/as 
done. Thus, the evidence indicates casual environmental factors liad 
no influence on the consistency or computational time of a solution. 

A. Testing With Pre-Determined Points. Accuracy checks were 
accomplished by picking ra^iidom locations in the Loran-C coverage 
region for a given triad and computing the corr<^*«,)onding TD readings 
for that triad of stations. These TDs were theu used as input to the 
coordinate converter and the results were compared with the original 
position coordinates. These are shown In Figures 5-1 and 5-2. Figure 
5-1 are the results obtained using an early version of the program 
which employed a spherical earth model, while Figure 5-2 shows the 
results of the latest version which uses an ellipsoidal earth model. 

In both cases, the center point used in the TD-to-position calcula- 
tions was chosen to be the center of the l^ran-C triad used. This is 
done by finding the average latitude and longitude of the master and 
the two selected secondary stations used. Comparing the two results 
shows that with t)ie spherical approxiixation, errors tend to increase 
sharply as the distance from the center point increases beyond 150 
nautical miles. In the elliptical case, the errors are much smaller 
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Figure 5-1 . Test of Microcomputer Coordinote Converter 
with Sphericoi Model. 
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Figure ^2. Continued. 
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out to the edge of the coverage region and the reeultii tend to be more 
stable as a function of dlst/mce fron the center of the coverage 
region. 


The results presented here are nuch the sane In terms of 
accuracy as the tests done in Chapter III using the FORTKAN-*IV test 
program. The results in Figures 3*1 and 3*2 also Include the range 
and bearing angle calculations presented in Chapter Tli.D. The range 
and bearing calculations allow only a single wa/poiut to be stored in 
memory^ which must be established before TD-to*posit Ion updates are 
run. There are 11 steps required for each rangn/bearing computation 
and this is normally done immediately after the lat/long calculations. 
The output of the microcomputer program consists of the computed 
lat/long and the range/bearing stored in memory in both the binary 
floating-point format and in BCD for data storage and display. 

Addillonal testing of the microcomputer Loran*C system 
described in Chapter IV mu performed during actual flights in a Piper 
Cherokee and Douglas DC-3. These tests were done mainly to demon- 
strate the fact that the microcomputer system could be used in flight 
envlronMnts to provide navigation information. This flight testing 
also allowed the output of the coordinate converter to be recorded 
over a long period of time In flight situations for checking the sta- 
bility of the microcomputer hardware and software system. Approxi- 
mately 30 hours of flight testing were done and it was demonstrated 
that the output of the coordinate converter » in the form of range and 
bearing to a waypoint » could be useful to a pilot. Figure 3-3 shows a 
plot of the recorded data from the coordinate converter on a particu- 
lar flight along with the recorded data from a commercial I/>ran-C 
receiver. From here It can be noted that there is an obvious bias in 
the microcomputer data on the east-west track to the north; this track 
coincides nearly %d.th the LOP from the Seneca Carolina Beach pair, 
which was used for one of the TDs. The other pair used was Seneca- 
Dana» in the 99600 Loran-C chain. This particular flight path 
illustrates an interesting characteristic; the flight path along the 
north-south direction coincides more nearly to the Seneca-Dan.i LOP, 
which is a stronger signal pair, and the oast-wst track coincides 
with the weaker LOP signal pair. As is evident, the north-south track 
is on course and cast-west track is slightly north of course. The 
comoierclal receiver flown incorporates signal path corrections, wliile 
the microcomputer system does not. It is obvious that signal path 
corrections are needed, although the inaccuracy without them amounts 
to onc-half to one nautical mile. 

The flight testing was done to assess the microcomputer- 
controlled navigation system's capabilities rather than evaluate abso- 
lute accuracy. Analysis of other flight tests similar to those shown 
in Figure 3-3 showed typical position uncertainties of one-half to one 
nautical mile. For general aviation use, this represents sufficient 
accuracy for enroute navigation and to place the pilot near navigation 
fixes at an airport for landing approaches. Much of the error noted 
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Fiaure 5>-3. Results of Test Flight of Coordinate Converter with Commerciol Receiver. 
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In flight touting could be treced to the receiver where the TUh are 
derived. Aeeumlng that the TDe could be derived correctly, the 
eveiieble accuracy would be on the order of 0. 1 to 0.1 nn. TV.e 
effecta of the derivation of the TDa In the receiver in a function of 
how the pulae ia oanpled. 

Flight leating alao ahowcd that the Loran-C algnala imny be 
received well at low altitudes, even in the hilly terrain where th«: 
teata were conducted. Fffecta on the r.rnuncl siu li .*i< hO ||/. imwor-linr 
interference nre fmich leas prononnrcH during fllrlu. 

K.rr«»rs as Function of Dl atancc . Obvlouaiy, the choice 

of the central ;r>lnt uaed to compute the conatanta Cl through C9 in 
the aphcrlcal cai^e, and Cl through C14 in the clLipaoidai case, has a 
definite bearing on how accurate the solution %rill be ct any given 
point. If it ia kiown that moat operationa will take place within a 
amall part of a given triad coverage region, the central point could 
be choacn there. In general, though, the accuracy will deteriorate 
rapidly aa the dlatance from the central point eatends beyond 200 to 
230 milea. The central point in the microcomputer program was choaen 
CO correapond to the center of the triad region to provide good 
accuracy throughout the triad’ a coverage region. 

Fo**ib Ie Fixe s for ^rorjj^. Certainly, more accurate models 
of the eartW proTrie over whic^i the Loran-C algnal travels would pro- 
vide more accurate aolutlona. This becoiacs a very difficult task in 
practice because of al? the various profiles which exist in any ‘;lven 
area and it ia not practical to store these and i^erfora the necessary 
geodesic calculations on a small microcomputer syotem. 

Many modern receiver designs are incorporating methods of 
applying propagation corrections to Che received time differences as 
the position calculations are made. This la done by going to various 
locations and observing actual time differences. Propagation correc- 
tions become very important in determining how accurate any receiver 
•yatam la, although aa mentioned in Chapter II, these correct iona typi- 
cally are on the order of several tenths of a mile, thus ''he 
complexity of including propagation correct ns may not be Justified 
for some situations. It is alno possible to include minor earth M>del 
corrections into the propagation corrections, although these would, in 
most cases, be very small corrections. 

Kb was discussed «bove relating to Figure 3-3, the best way 
to provide more accurate TD-to-poiltlon calculations is to incorporate 
soma sort of method of compensatl.\g for weak signals from the Loran-C 
•tations. The strength of the signal Is generally related to the 
signal path length and can thus provide useful information for other 
corrections, such as the profile of the earth over which the signal 
travailed. Although this can decrease inaccuracies from one nautical 
mils to less than onc-tanth nautical mile, the added complexity of 
getting station aignal-lo-noise ratios from the Loran-C receiver and 
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applying corrections basod on these can sometiaes be too luuch of a 
burden in processing tlae or storage requlreaentu. It may be quite 
feasible to have several ground monitoring stations in the coverage 
region to determine local corrections (differential Loran) which could 
then be disseminated by some means to the pilot or operator and 
entered into the Loran-C navigation system* 
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VI. CONCLUSIONS AND RECOMHEN DAT IONS OF POOR QUALITY 

In thl8 paper* one objective has been to show the suitability 
of using Loran-C as a navigation systea for general aviation use. 
Loran*C has many inherent capabilities which aakes it attrar.tive with 
respect to along and cross-track accuracy and signal coverage. Key 
conclusions reached as a result of this research are: 

1. The obstacles to fully implementing a complete Loran-C navi- 
gator because of the hardware and the computational burdens are demon- 
strated to be removed with the wide availability of lo%#-cost 
microprocessors and related devices. 

2. The problem of converting Loran-C data to a more natural 
and easily understood form for the pilot using simple and fast tech- 
niques has been solved. 

3. The cost of implementing a microprocessor-based Loran-C 
navigator unit can be comparable to the cost of a small microprocessor- 
based* general-purpose computer system in terms of hardware and soft- 
ware required. 

4. The TD-to-position algorithm developed has comparable 
accuracy to that which can be achieved with raw Loran-C data; approxi- 
mately 0.1 to 0.3 nm or less. 

This work effort lias concentrated on the computational aspects 
of a Loran-C navigator. The fact that transcendental functions are 
involved has impeded Implementation ot these ideas in the iiast; these 
functions are still difficult to implement on an eight-bit micropro- 
cessor. The tradeoffs between extra liardware to handle these func- 
tions or doing them in software must be made in terms of speed since 
other tasks must be done in the navigator. Thus* use of an external 
mathematics chip* such as the Am 9311* has a strong advantage in speed 
over software mathematics routines. This is especially true in light 
of recent price reductions for such satellite processors. 

Pecause of velocities typical of aircraft operations* it Is 
ideal to have complete position updates in one second. It should be 
noted that extra area-navigation computations in addition to the 
Loran-C position calculations (stated at 0.18 second) may exceed this 
desired value. The Am 9511 may be run asynchronously at speeds up to 
2 MHz allowing the microprocessor cime for other navigator tasks. 

In view of the speed and accuracy of processor loading parame- 
ters of the coordinate conversion routine* the conclusion is that 
these techniques may be successfully used for general aviation Loran-C 
navigation application. 

During the course of the investigation for this paper* several 
i ems were found which need to be addressed further. One is the 
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acceptance of U>ran*C as a primary means of navigation for aircraft. 
The present Loran-C system is maintained primarily for the maritime 
interests which means that some engineering changes may be necessary 
before full implementation of Loran-C for aviation could take place 
[37,38] • additional stations are desired to cover fully the United 
States, and it is important to place these to provide maximum capabi- 
lity to users, especially those in areas of rough terrain or areas 
lacking in numbers of other services [39,40]. 

A very important issue is the consideration of how over-land 
propagation of the Loran-C signal affects the phase velocity of the 
signal. Explicit solutions to this are difficult to implement because 
of the very inconsistent nature of the propagation path. Empirical 
models are usually tsed. This is done by isonitoring the Loran-C 
signal at many locations and averaging the results over long periods 
of time. A table of results could be published and the values 

could be stored iu the navigator unit, or broadcast to the aircraft 
much r.s barometric pressure readings are. Another Important effect on 
the propagation path model is the g«^odctic model of the earth [41]. 

It is felt that these problems are manageable [42] and are not a 
serious impediment to fully using Lova^'-C as a primary means of air 
navigation. 


ORIGINAL 

POOR QUALITY 


- 80 - 


Vll. 


acknowleix;kments 


Acknowledgeaeat is made Co the following people %dio aided with 
Che research reporced in this paper: Dr. Roberc W. Li 1 Ley, assoc iaCe 

direccor, Hr. Ralph W. Burhans, projecC engineer, and sCudenC 
researchers Daryl HcCall, Sceven Yosc and Scanley Novacki. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


- 81 - 


ORIGJMAL PAGE IS 

Vm. REFERENCES OP QUALITY 

[1] . Kaycon, Myron and Walter R. Fried, "Avionics Navigation 

Syateas," John Wiley and Sons, Inc., New York, NY, 1969, p. 160. 

[2] GPS Special Issue, Navigation, Vol. 25, No. 2, Summer 1978. 

[3] . Pursel, R. H. , "Flight Test Investigation of LORAN-C for Civil 

Aviation Applications," Navigation, Vol. 27, No. 3, Fall 1980. 

[4] . HcConkey, E. D. , "Evaluation of LORAN-C for Non-Precision 

Approach Applications," Navigation, Vol. 27, No. 3, Fall 1980. 

[5] . Natarajan, Krishnan, "Testing of Loran-C For General Aviation 

Aircraft," NASA Conference Publication 2176, Proceedings of 
Joint University Program For Air Transportation Research - 1980 
Conference, December 11-12, 1980. 

[6] . Wong, Gene A., "Analysis of Loran-C System Reliability For 

Civil Aviation," Proceedings of Thirty-Seventh Annual Meeting, 

The Institute of Navigation, Annapolis, Maryland, June 1981. 

[7] . Stremler, Ferrel C. , "Introduction to Communication Systems," 

Addison-Wesley Publishing Compary, Reading, Mass., 1979, p. 509. 

[8] . Radio Amateur's Handbook, 5lai. edition. The American Radio Relay 

League, Newington, CT. , 1974, pp. 559-66. 

[9] . Swanson, E. R. and F. C. Roble, "VLF Navigation," NELC Technical 

Document 233, Naval Electronics Laboratory Center, San Diego, 
Cal., January 1973. 

[10] . Lepanto, Janet A., "Omega Coordinate Transformations and 

Propagation Corrections," MS Thesis, Princeton University, 
Princeton, N. J. , September 1980. 

[11] . Chamberlin, Kent A., "Digital Correlation Detector for Low-Cost 

Omega Navigation," M.S. Thesis, Ohio University, Athens, Ohio, 
December 1975. 

[12] . Pierce, J. A., and R. H. Woodward, "The Development of Long-r«ange 

Hyperbolic <^vigation in the United States," Navigation, Vol. 18, 
No. 1, Spring 1971. 

[13] . Hefley, Gifford, "The Development of Loran-C Navigation and 

Timing," NBS Monograph 129, U. S. Dept, of Commerce, National 
Bureau of Standards, Boulder, Colorado, October 1972. 

[14] . Snerman, H. T. and V. L. Johnaou. 'The Loran-C Ground Station," 

Navigation, Vol. 23, No. 4, Winter 1976-77. 

[15] . "Loran-C User Handbook," COMDINST MI6562.3, Department of 

Transportation, United States Coas‘ 'luard. May 1980. 


- 82 - 







[ 


i 

i' 

1 





[16] . “CoordlnACet and Ealsaion Delays of Loran-C Stations,” Tine 

Service Announceaent, Series 9, No. 150, U. S. Naval Observatory, 
Werhington, DC, Septeaber 1979. 

[17] . Op clt., Kayton and Fried, pp. 26-7. 

[18] . "Method of Coaputlng Loran Hyperbolic Lattice,” Defense Mapping 

Agency, Hydrographic Center, Washington, D. C. , 1977. 

[19] . Uolastrom, J. S. , "A New Approach to the Theory of Geodesics on 

an Ellipsoid,” Navigation, Vol. 23, No. 3, Fall 1976. 

[20] . Saaaddar, S. N. , "The Theory of Loran-C Ground Wave Propagation - 

A Review,” Navigation, Vol. 26, No. 3, 1979. 

[21] . Martin, 0. L. , Private Coonunicatlon, Defense Mapping Agency, 

Washington, D. C. , Septeaber, 1979. 

[22] . Op clt., Kayton and Fried, pp. 27-8. 

[23] . Caapbell, A. C. , "Loran to Geographic Conversion and Geographic 

to Loran Conversion,” Informal Report N-3-64, Naval Oceanographic 
Office, Washington, D.C. , June 1964. 

I 

[24] . Piecuch, L. M. and R. W. Lilley, "Interactive Loran-C to 

Geographic and Geographic-to-Loran-C Coaputation,” NASA TH 52, 
Avionics Engineering Center, Ohio University, Athens, Ohio, 

August 1977. 

[25] . Friedland, B. and H. F. Hutton, "New Algorithms For Converting 

Loran Tiae Differences to Position,” Navigation, Vol. 20, No. 2, 
Summer 1973. 

[26] . Roth, W. F. , A. R. VanKoevering and D. Wozniak., "Hyperbolic 

Coordinate Converter,” Lear Slegler Publication No. GRR-65-1023, 
Lear Slegler, Inc., Grand Rapids, Michigan, January 1965, pp. 3-1 
to ^6. 

[27] . Rasin, Sheldon, "Explicit (Noniterative) Loran Solution,” 

Navigation, Vol. 14, No. 3, Fall 1967. 

[28] . Fell, Harriet, "Comments on Loran Conversion Algorithms," 

Navigation, Vol. 22, No. 2, Summer 1975. 

[29] . Op cit. , Roth, W. F. , et. al., pp. 3-7. 

[30] . HCS6500 Microcomputer Family Hardware Manual, Publication Number 

6500-50, MDS Technology, Inc., Norristown, PA, January 1976. 

[31.1. MCS6500 Microcomputer Family Programming Marj&l, Publication 

No. 6500-50, MOS Technology, Inc., Norristown, PA, August 1975. 

[32]. Lilley, Robert Vi., "An Assembler For the MOS Technology 6502 
Microprocessor as Implemented in Jolt(TM) and KIN-l(TM)," 

NASA IM 44, Avionics Engineering Center, Ohio University, Athen:i, 
Ohio, November 1976. 


-83- 




i 


1 




i 


w 


[33] * Llllay, R. W. and D» L. McCallt "A Lorao-c Prototype Navigation 

Racalvar For General Aviation,” OU/NASA TH 80, Avionics 
Eaginaerlng Center, Ohio University, Athens, Ohio, August 1981. 

[34] . Fischer, J. P. and R. U. Lllley, ”A Video Display Interface For 

Che Loran-'C Navigation Receiver Developaent fasten, ” OU/NASA 
TM 38, Avionics Engineering Center, Athens, Ohio, Hay 1978. 

[33]. "N6800 Hicrocoaputer Syatea Design Data,” Motorola Semiconductor 

Products, Inc., Phoenix, Arizona, 1976. 

[36] . Bucklen, Willard, John Eldon, Louis Schlrm and Fred Wllllans, 

”Slngle-Chlp Digital Ihiltlpllcrs Fora Basic DSP Building Blocks,” 
EDN, Vol. 26, No. 7, April 1, 1981. 

[37] . Federal Redlonavlgatlcn Plan, Volume 1, Radlonavlgatlon Plans 

and Policy and Volume 2, Requirements, Departments of Defense and 
Tranaportation, Washington, D. C. , July, 1980. 

[38] . Tyler, J. S. , and F. C. Karkallk, ”Redlo Navigation Systems - 

Current Status,” North Atlantic Treaty Organization Advisory 
Croup For Aerospace Research and DevelopMnt, ACARlV'AC-243 
(Principles and Operational Aspects of Precision Position 
Determination Systems), NeulllySur-Seine, France, July 1979. 

[39] . Polhemus, William L. , ”Operational and Economic Benefits Deriving 

From Use of Loran-C RNAV,” Proceedings of the Thlrty-Se— nth 
Annuzl Meeting, The Institute of Navigation, Annapolis, Maryland, 
June 1981. 

[40] . Rzonca, Lorraine, ”Loran-C RNAV in Mountainous Areas,” 

Proceedings of the Thirty-Seventh Annual Meeting, The Institute 
of Navigation, Annapolis, Maryland, June 1981. 

[41] . Kuebler, W. and S. Sommers, ”The Role of Che Earth's Shape In 

Navigation: An Example/* Navigation, Vol. 28, No. 1, Spring 1981. 

[42] Walker, U.L. and R. Ellerbe, ”Deslgn Study Report for General 
Aviation Loran-C receiver,” Teledyne Systems Company, U.S. 
Department of Transportation, Federal Aviation Administration, 

DOT FAA/RD-81/36, Washington, D.C., May 1981. 

[43] . IBM Syatem/360 and System/370 FORTRAN-IV Unguage, IBM 

Publication CC28-6513-10, May 1974. 


original page is 
OF POOR QUALITY 


- 84 - 


IX. 


Ar‘a.rOICES 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Appendix I . Program Hating for position-to-TD 
prediction program. This program is written in standard 
FORTRAN-IV computer language for the IBM/370-158 [431. 
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- THIS PHOGRAM IS A TEST PHOGRAM FOR CCMPUTING L0RA>«-C 

• TIME DIFFERENCES. THE EQIMTION USED IS: 

- TO • BETA ♦ DELTA ♦ TS - IN 

• WHERE BETA IS THE BASELINE LENGTH, AND DELTA IS THE 

- SECONDARY COOING DELAY. TS AND TM NUST BE CALCUUTED 

- AND REPRESENT THE BASELINE LENGTH FR(M RECEIVER TO THE 

- SEOOIOARY STATION AND fHE BASELINE LENGTH FROM T>€ 

- RECEIVER TO THE MASTER, RESPECTIVELY. THESE TNO ARE 

- CCMPUTED BY "ARC* AND THE ACCURACY DEPENDS ON THAT 

- SUBPROGRAM. ALL GEODESIC CONSTANTS ARE NM) 27. 
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C - J.P.FISCHER 07,79 

C - 

C - CHAIN CONSTANTS ARE FOR U.S. EAST COAST (9960). 

C - 

REAL'S PHIR,GAMR.PHIM,GAMM,PHIW,GAMN.PHIX,GAMX,PHIY,0(MY,PHIZ.GAMZ 
COMMON PHIR.GAMR 

DATA PHIM,GAMV0.749S002761.1.34('«70724/ 

DATA PH IN ,6AMN/0.81 69491 990.1.1 89999303/ 

DATA PHIX,6AKX/0.7200063971, 1.221349098/ 

DATA PHIY,GAMir/0.9949097338, 1.399840319/ 

DATA PH IZ.GAMZ/0.6999476439, 1.926928009/ 

DATA C06N/ 1 . 379724E4 / .CD8X/2. 696991 E4/ , 008 Y/4 . 222 16 1 E4 / ,C06Z/S .716 
>209E4/ 

C - 

9 PRINT 21 

REA0(9, 1 .EN0>6) I0E6,MIN.SEC, I0E61.MI N1 ,SEC1 
PH IR-( I KG*i Ml N»SEC/60. )/60. ) *3. 14 1 99269/18C. 

GA»«< 1 0EGHCMI NI^SECI/GO. )/60. ) *3. 141 99269/180. 

TDM-AFC(PH IM.GAMM) 

TON<DeN^ARC(PHIW ,GAMW)-TDM 
TOK-COeX^ARC(PH I X.GAMX >-T[M 
TDY-C06YMHC(PH I Y,GAMV)-T1>1 
TOZ -COB Z4^ARC ( PH I Z ,GAMZ )-T0M 
PRINT 22,TON,TOX,TOY,TDZ 
GOTO 9 
6 STOP 

1 F0RMAT(I4,IX,I2,IX,F9.0) 

2! F0f«4AT('0ENTER RECEIVER LATITUOE/LONGITUOE. V S(XX) MM SS.SS* ) 

22 F0»»4AT(' N - ',™.2,9X,'X - •,F8.2,SX,»Y • ',F8.2,9X,'Z - ',F8.2) 
END 

FUNCTION ARC(PHII,XLNG2) 

C 

C THIS PR0GRM4 IS DERIVED FROM KAYTON AND FRIED; 

C IT IS USED TO CALCULATE THE BASELINE DISTANCE BETWEEN TWO 

C POINTS ON THE REFERENCE ELLIPSOID BY DER|Vr« THE GEODESIC 

C ARC LENGTH IN A PLANE TH10UGH THE C^EMTER OF THE ELLIPSOID. 

C REFERENCE ELLIPSOID USED IS CLARKE 1966 (NAD 27). 

C 

C JP FISCHER 7/79 
C 

implicit REAL*8(A-H,0-Z) 

realh arc 

COMMON PHI.XLNGI 

DATA A/6,378206406/,F/3.3900793040-3/ 

DATA a/299. 7924 98/ 

C 

OXLNO«OABS(XLNC I-XLNG2 ) 

BETAOATAN( (I .-F)*DTAN(PHI ) ) 

BETA I K)ATAN( ( 1 . -f ) *DTAN(PH 1 1 ) ) 

CI-OOOS(8ETAI )«DSI N(OXLNG) 

C2-0C0S(BCTA)*DSI N(BETAI )-DSIN(BETA)'DCOS(HETAI )*DCOS(DXLNG) 

C3-0SI N(BETA)»DSIN(BETAI )♦OCOS(BETA)•OCOS(BETAI )»OCOS(OX1 MG) 
PSI-0ATMKC1/C2) 

THETA»OATAN( (C2*ttOS(PS I )*CI *OSIN(PS I ) )/C3) 

XM>(OSIN(BETA)-FDS INIBETAI ) )**2 
XN-( (OS IN(BETA)-OS IN(BETAI ) )/DSI N(THETA) )«»2 
XU»(1.-0C0S(TMETA) )/DS INITHETA)'(THETA-0SIN(THETA) )/DS IN( THETA) 
XV»( 1 .•K)COS(ThETA) )»(THCTAM)SI N(THTTA) ) 

RHCK>A8S( A'TH£TA-A«F« (XM'XU^XN*XV)/4. ) 

ARORHCVa 

return 

CIO 
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Appendix II» Progran listing for TD-to-lat/long 
conversion program. This program use written In 
standard FORTRAN-IV programming language and run on 
the IBM/370-158. 



C THIS PRQQRAM USES AN EVLICIT NON-ITEKITIVE APPROACH TO SOLVE THE 

C LORAN SOLUTION OF CONVERTING TINE DIFFERENCE REM)INCS TO LATITUDE 

C LONGITUDE. 

C EARTN CONSTMTS NAO-27 (CUfKE 1866 SPHERICAL EARTH MDOEL) 

C 11/1980 J.P. FISCHER 

C ELLIPSlOO CORRECTIONS APPLIED 6/1981 

C 

DATA TCT,TC2/3.9E-2,5.4E-2/ 

DATA T»MY/0.1SI29298/,T»«(Z/0.U6499S7/.»«</46.986746/ 

DATA C7Wr/0.988S7709/.S1Mr/0. l907l607/,C1MZ/0. 98899478/ 

DATA S1MZ/0.l479M43/.Q(K/0.37466568/,SXK/0.927l6O79/ 

DATA CI/2.983207I/.C2/4.311 1689/ ,C3/*/. 17171 16/ 

DATA C4/4.0896360/.C9/-4.4643647/.C6/1. 1496479/ 

DATA C7/.9.0910869/.C8/9. 1 189063/.C9/.9992 1480/ 

DATA CIO/1. 9830077E>3/ .Cl l/2.9836473E-4/,CI2/-2.96r03l8£-3/ 

DATA Cl 3/1. 3890077E*3/,C14/1 .0034014/ 

C 

7 PRINT 12 

READ<9,3,ENI>I0)TY,T2 
TY»TY*l,E-6 
TZ«TZ»1.E-6 
C 

PY^XNR* (1 Y-TCY 1-TI9IY 
PZ-XNR* (T2-TCZ) 'Tl*« 

CPY<OS(PYl 
SPY-SI N(PY) 

CPZ¥X)S(PZ) 

SPZ-SIN(PZ) 

AY-(CPY-CTNV)/STNY 
AZ-(CPZ-CTNZ)/STNZ 
BY-SPY/STMY 
BZ-SPZ/STMZ 
U1-AY»CXK-AZ 
U2-AY«SXIC 
U3-AZ*BY-AY»BZ 
Ul^U1»U1^U2nj2 

C08Y-(U3*U1'HJ2«SORT(UU-03*U3 ) )/0U 
Tl4«S-ATAN(AY/(BY«C0eY) 1 
CB<0S(TtMS) 

CA<0S(T»«*SVPY) 

COCOS! TlftS-iPZ) 

c 

F-C1*CA^«C8>C3«CC 
OC4»CAK9»C8406»CC 

h-C7*cah:8*>»r: 9» cc 
c 

TH6S-ArAN(«M:iO)/IF4C1 I) ) 

P»CS-ATAN(CI 4 »C1 4*S I N( THGS ) • (H«CI 2) / (G«^C 1 3) ) 

CALL LLR0(PHBS,THGS) 

GOTO 7 
10 STOP 

3 FORMATIFIO.O) 

12 FORMAT! ENTEH TIME-OIFFERENCES.' ) 

EM) 

SUBROUTINE ROLL!PHI .THE,* ) 

C 

C THIS SUBROUTINE CONVERTS GEOCENTRIC COORDINATES ENTERED BY THE 

C USER TO RADIAN COORDINATES. IM>UT FORM IS: DODO MM SS.SS 

C NHERE '0000' IS THE DEGREES PORTION OF THE LAT. OR LONG., IICLUDIN 

C SIGN. MM IS THE MINUTES PORTION, AND SS.SS IS T>€ SECONDS PORTION. 

C READ FOlWAT IS: M, IX, 12. tX.Fj.O. 

C 

implicit REAL*8!A-H,0-Z) 

real 44 PHI, THE 

DATA PI/3.1413926939898/ 

DATA MSC1/'LATI'/,MSG2/*TU0E'/ .MSG3/': '/.MSG4/' L0*r.*/,MSG5/' I TU 

>0'/,*«06/'E: •/ 

PH-PI/180. 

c 

c PRorr USER, 
c 

PRINTI,MSCI.MSG2.MSC3 


OWOJNAL ^ 

OF POOR QUALITY 
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REAO(3,IO,ENO-2l>ID1.IMI,SSI 

PHI«SMr^(PI1«(0fL0AT(l0t)*(0fLQAT(l»4l HSSI/60. )/60. ) ) 

PROMPT USER FOR LONGITUDE ENTRY, ORIGINAL PAGE i3 

PRINTI,MSG4.MSG»,MSG6 POOR QUALIT'^ 

READ(^,l0,END>21)l01,IM1,SSt 

TNE-SNGL(PI I • (0FLQ^T{ 1 01 )>(0FLQAT( IM1 HSS1/60. )/60. ) ) 

return 

21 RETIFN I 

ENTRY LU®(PHI,TME) 

THIS SUOSCCTION DOES THE REVERSE OF THE ABOVE GIVEN RADIAN 
CCOROINATES. 

PI2-180./PI 

PHI^I2*PHI 

THEi^l2nH£ 

LDt«PMI 

»4I"(PHI-L01)«60. 

LM1-XM1 

551- (XMI>LMI)*60. 

L02-THE 

m2-(THC>L02)*60. 

W>XM2 

552- (XM2-LM2)*60. 

PRINT3,L0M>«1,SS1,L02,LM2,SS2 

RETURN 

1 FORMAT! ' ENTER ' ,3A4/» 0000 MM SS.SS') 

10 F0I«4AT(I4,1X, I2,1X,F9.0) 

5 FORMAT!' LATITUDE ■ • ,14 , 1 X, 12, IX.F5.2/' LONGITUDE ■ ',14, IX, 12,1 

>X,F5.2) 

EW> 
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S?'pcOR 


Appendix III . Prograa listing for olcroprocessor 
version of TlHto-l^.t/long program. This program was 
«frltten In standard HOS6502 assembly language and 
assembled by a cross assembler on the IBM/ 370-150. 



THIS PWOORAM IS oeSIONEO TO CONVERT LORAM-C T IME-0 IFFEREHCFS 
TO LATITUOC/LONGITUOE AM) RANGE/BEARI NO TO A MAYPOINT 
ON THE MICKCOOHPUTER USING THE AM95HA MATH CHIP. THE 
ALG^MTHM IS BASED ON FORTRAN PROOR/M •DCXLRN,* BCD-1 >HEX 
AND HEX-TO-8GD CONVERSIONS ARE M/V>E AT THE BEGINNING AND 
END OF PROGW EXECUTION; ALL INTERNAL CALCULATIONS ARE 
MADE USING BINARY FLOATI NG-PO INT. ^L SUBROUTINES ARE AT 
THE END OF THE MAIN PROGRAM. THE N«J4BER TABLE AREA IS 
DESIGNED TO BE PLACED AFTER THE SUBROUTINES; CONSTANTS ARE 
FIRST, CALCULATED VARIABLES LAST. 

VERSION 1.0, SPHERICAL MODEL, 11/1980, J.P.FISCMER 
VERSION 2.0, ELLIPSOID MODEL, 7/1961, J. P. FISCHER 


PIAA EQU $9000 PERIPHERAL AND DOR SIDE A 

PIAB EQU $9002 PERIPHERAL AND DOR Sl^iE <3 

• 

* /949S1IA COMMANDS. 

• 

FAOD EQU $10 

FSUB FQU $1 1 

FMUL EQU $i2 

FDIV EQU $13 

SORT EQU 1 

SIN EQU 2 

COS EQU 3 

ATAN EQU 7 

PTOF EQU $17 

PUPI EQU $1A 

FLTD EQU $1C 

FIXD EQU $1E 

• 

* VARIABLE NAME TABLE 


• THE FOLLONING ARE CONSTANTS. 

• 

TCY EQU 0 

TCZ EQU TCYM 

THMY EQU TC244 

jmz EQU 

xm EQU T^9«Z44 

CTMY EQU XNR4A 

STMY EQU CTNY-M 

CTMZ EQU STMY^4 

STMZ EQU CIMZM 

CXK EQU GTHZ^M 

SXK EQU CXK44 

Cl EQU S)CKM 

C2 EQU CIM 

C3 EQU C2M 

C4 EQU C3^4 

C3 EQU C4M 

C6 EQU CV4 

C7 EQU CtH 

C% EQU C744 

C9 EQU C6M 

CIO EQU C9M 

Cl 1 EQU C10M 

C12 EQU CUM 

C14 EQU C12^4 

EM6 EQU C14M 

RCR EQU EM6M 

TNO EQU RCIH4 


ORIGINAL PAGl TJ 
OF POOR QUALITY 
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C256 

EQU TWM 


PI80 

EQU C25644 


C60 

EQU PiaOM 


PHI 

EQU CtOH 


RLAM 

• 

EQU PHI44 


• 

• 

THE FOLLOW INC ARE CALCULATED VARIABLES. 

TY 

EQU RLAMM 


TZ 

EQU TY^4 


PY 

EQU TZ44 


PZ 

CPY 

EQU PY^ 
EQU PZ^4 

OWGINAL PAGE O 

SPY 

CPZ 

EQU CPYM 
EQU SPY44 

OF POOR Q'-^’ 

SPZ 

EQU CPZ 44 


AY 

EQU SPZ44 


AZ 

EQU AYM 


BY 

EQU AZ44 


6Z 

EQU BY44 


Ut 

EQU BZ^4 


U2 

EQU Um 


U3 

EQU U244 


UU 

EQU U>M 


C08Y 

EQU UUM 


jms 

EQU C08Y 4 


CB 

EQU T)94SM 


CA 

EQU C844 


CC 

EQU CA44 


F 

EQU CCM 


G 

EQU F44 


H 

EQU Gf4 


THGS 

EQU 


PHGS 

EQU THBS44 


CFl 

EQU PHGS^4 


DLAM 

EQU CFH4 


OPHI 

EQU 0LAIH4 


TEMP 

• 

EQU DPNI44 


• 

• 

PAGE-ZERO ASSIGNMENTS 


ORG 0 


TCM 

BSS 4 

PACKED BCD /W TENTH DIGIT 


BSS 4 
ORG SBE 

PACKED BCD /W TENTH DIGIT 

AXF 

BSS 1 


AOOB 

BSS 1 


6ASEI 

BSS 2 

BASE AOCY^ESS FOR MMBER MOVE 

XLIM 

BSS 1 

INDEX LIMIT USED IN NUMBER CONVERSION ROUTINES 

XJDfP 

BSS 1 

SAVE AREA FOR X-REGISTER 

YTEMP 

BSS 1 

SAVE AREA FOR Y-REGISTER 

COUNT 

BSS 1 

COUNT-DOWN REGISTER 

VY 

BSS 1 


flag 

BSS 1 


Ht^ 

BSS 3 

table for BUILDING up hex NlfMBER FROM BCD 

RES 

BSS 3 

RESIDUE TABLE FOR BINARY MULTIPLICATION 

OVDN 

OSS 2 

HEX NUMBER TO BF CONVERTED TO BCD 

ORES 

nss 5 

RESULT OF HEX-TO-BCD CONVERSION 

OYSR 

BSS 1 

CONSTANT USED IN BINARY DIVISION 

RMNGR 

BSS 1 

REMAIIOER FROM BirUUTV DIVISION 


The following are the 

QOMPUTED LaTITUDE/LONGITUDE AM) RANGE/BEARING. ALL 
ARE IN PACKED BCD fOmhJ. T>€SE ARE DESIGNED TO 
INTERFACE TO SENSOR SOFTWARE. 


LAT 

RHO 


BSS 3 
BSS 3 
BSS 2 


FORMAT: (EX.) 39 19 20 - DEGREES, MINUTES, SEOON 
FORMAT: (EX.) 02 05 56 - DEGREES, MJNUTES, SECON 
FORMAT: (EX.) 00 96 - NAUTICAL MILES /W TENTH 01 


ess 2 FORMAT: (EX.J 22 04 - OEUREES /K TENTH DIGIT 

BASE 83 S 2 ADDRESS OF FLOATING-POINT TABLE 

OhiuiiiAL PAar F3 

OP POOR QUALiry 


• BGD-TO-HEX CONVERSION 

• 

TOMDC LOY -6 
MM -0 
STA HEX 
STA MEX^I 
STA HEX>2 
STA FLAG 
J»f> LFOiF 
UFOif) OeC FLAG 
LOA TQA,X 
LSR A 
LSR A 
LSR A 
LSR A 
JMP T05 

LFQVF \fC FLAG 
LOA TQA,X 
AND -SF 

I fitt 

TQ5 CLC 

AOC HO>2 
STA HEX4«2 
BCC NX 
lie HEX^I 
NX OEY 

BEQ TX 
STY XT&P 
m 

• kjltiply partial sum by ten 

• 

LOA -0 

STA RES CLEAR MULT. TABLE 

STA RES^I 
STA RES42 

LOA -10 DIVISOR 

LOY -8 
T02 CLC 

ROL RES^2 
ROL RES^l 
RX RES 
ASL A 
BX N0C2 
PKA 
CLC 

LOA RES-^2 
ADC HEX>2 
STA RES^2 
LOA RES^I 
AX I€X4>1 
STA RES^I 
LOA RES 
AX HEX 
STA RES 


SIX DIGITS 

00 LONER PART FIRST 
GET A BYTE 

X THE MAIN INVERSION 
GET A BYTE 

REMOVE UPPER FOUR BITS 

FOR ADDITION 

ADO DIGIT TO PARTIAL SUM 
GO IF NO CARRY OUT 

NEXT DIGIT 
IF DONE, LEAVE 
SAVE COUNTER 


ORG II7X 

XNYERT TO'S FROM BX TO FLOATING-POINT 
LOX -0 

LOY *TY XNVERT TOA FIRST 

STY YTEMP 
LOA -2 

STA COUNT CONVERT TNG TOS 

CLD SET DECIMAL MODE OFF 





i 



N0C2 




D€Y 

BN£ TQ2 
LDY xreNP 
LOk RES 
$TA HEX 
LOA RES^t 
STA HEX^I 
LOA RES42 
STA HEX^2 
LQA FLAG 
BNE UFOUR 
JMP LFOOR 


ORIGINAL PAGE IS 
OF POOR QUALITY 


NON CHANGC IMTEG£R PART TO FLOATING-POINT 
THEN M)0 IN THE FRACTIONAL PART. 


TOO 


SUTO 






LOY YTENP 
STX XTEMP 
LOA -0 

STA (BASE),Y 
IMY 

LQA HEX 

STA (BASE),Y 

INY 

LOA HEX^t 
STA (BASE),Y 
INY 

LOA HEX^2 
STA (BASE),Y 
LOY YTEMP 
JSR PUSH 
LOA -FLTD 
JSR CMNO 
LOY -D46 
JSR PUSH 
L1X^ -FHJL 
JSR cmo 
LOY YTEMP 
JSR POP 
DEC COUNT 
BEQ SUTO 
LQK XTEMP 
INX 

LDY «TZ 
STY YTEMP 
JMP TOMEX 


SAVE THE CURRENT DIGIT LXATION 

CLEAR THE UPPER VHO OYTtS 

MSB OF HEX INTEGER 
PUT I N TABLE FOR 95 II 

LSB Of HEX INTEGER 


POINT TO TO NIXBER 
GIVE IT TO 951 1 

CONVERT INTEGER TO FLOATING-PO INT 
CONSTANT IE-7 

CONVERT FROM MICROSECONDS... 

TO SECONDS 

GET TIC TD LXATION AGAIN 

AND STORE THE TD 

SEE IF BOTH TDS CONVERTED 

IF SO, START TD-TO-POSITION CONVERSION 

SET UP FOR NEXT TO 

ADDRESS OF TDB 

store it 

REPEAT 


CALCULATE 'PY* 


LDY -XNR 

•XNR* 

JSR PUSH 

“'SH ON STACK 

LOY -TY 

•TY‘ 

JSR PUSH 
LOY -TCY 

•TCY* 

JSR PUSH 
LQA -FSUB 
JSR CMNO 

SUBTRACT TY-TYC 

lGA =FHiL 
JSR CMNO 

XNR«(TY-TYC) 

LDY mTrny 

•TIMY* 

JSR PUSH 
LQA -FSUB 
JSR CMNO 

XNR«(TY-TZ)-Tlf4Y 

LOY M»Y 

•PY* 

JSR POP 

PUT *PY* INTO TABLE 


CALCULATE *PZ • 


LDY -XNR 'Xlfli 


C0002080 
C0002090 
C0002100 
000021 10 


C0002770 

CO0027B0 

C0002790 
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JSR PUSH 

PUSH OMTO STACK 

LDY -TZ 

•TZ' 

JSR PUSH 

PUSH ONTO STAa 

LOY -TCZ 

•Ta» 

JSR P'JSH 
LOA -KSUB 

PUSH ONTO STACK 

JSR OM) 
LOA 

TZ-TCZ 

JSR cmc 

XNR*(T2-TCZ) 

LOY mjmz 
JSR PUSH 
LOA -Fsue 

•Ttf»2' 

JSR CMNO 

XNR•{T^-TCZ)-T^♦^Z 

LOY -PZ 

LOCATION FOR 'PZ* 

JSR POP 

STORE RESULT IN 'PZ 

CALCULATE CPY,SPY,CP2,SPZ 

LOY ^Y 

•PY» 

JSR PUSH 

PUSH IT 

LOY -PY 

•PYi 

JSR PUSH 
LOA -COS 

DUPLICATE STACK 

JSR CMNO 

COS(PY) 

LOY -CPY 

FOR CPY 

JSR POP 
LOA -SIN 

GET IT 

JSR CMNO 
LOY -SPY 

SIN(PY) 

JSR POP 

GET SPY 

LOY -PZ 
JSR PUSH 

•PZ» 

LOY 

JSR PUSH 
LDA -COS 

»PZ* 

JSR cmo 

COS(PZ) 

LOY -CPZ 
JSR POP 
LOA -SIN 

LOCATION FOR CPZ 

JSR CMNO 

SIN(PZ) 

LOY -SPZ 
JSR POP 

LXATION FOR 'SPZ* 

CACULATE ' 

AY' 

LOY -CPY 
JSR PUSH 

'CPY* 

LOY -CTMY 
JSR PUSH 
LOA -FSUB 

•CTMY* 

JSR CMNO 

CPY-CTMY 

LOY -SIMY 
JSR PUSH 
LOA -FOIV 

*STMY* 

JSR O#© 

( CPY-CTMY) /STMY 

LOY bay 
JSR POP 

*AY* 

CALCULATE 

*AZ* 

LOY -CPZ 
JSR PUSH 

*CPZ* 

LOY -CTMZ 
JSR PUSH 
LDA -FSUB 

•CTMZ* 

JSR CMNO 

CPZ-CTMZ 

LDY -STMZ 
JSR PUSH 
LOA -FOIV 

•STMZ* 


ORIGINAL PAGE !S 
0? POOR QUALITY 
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JSR OMO 
LOY -AZ 
JSR POP 


(CPZ-CT>(Z)/S1WZ 
LXATION FOR *A^ 










t 


CALCULATE •BY* 


LDY -SPY 

•SPY* 

JSR PUSH 


LDY -SYMY 

•ST>IY* 

JSR PUSH 


LOM -POIV 


JSR CMN' 

SPY/STMY 

LOY -BY 

LXATION FOR •BY* 

JSR POP 

GET IT 

CALCULATE 

•BZ» 

LOY -SPZ 

•SPZ* 

JSR PUSH 


LOY -STRZ 

•STRZ* 

JSR PUSH 


LOA -roiv 


JSR CMNO 

SPZ/STXZ 

LOY ^iZ 

LXATION FOR •BZ* 

JSR POP 

GET IT 

CALCULATE 

•Ul* 

LOY »AY 

•AY* 

JSR PUSH 


LOY -CXK 

•CXK* 

JSR PUSH 


LOA -FRUL 


JSR CM40 

AY •CXK 

LOY -AZ 

•AZ* 

JSR PUSH 


LOA -FSUB 


JSR C»«0 

AY»CXK-AZ 

LOY 

LOCATION FOR *U1* 

JSR POP 


aACULATE 

•U2* 

LDY -AY 

•AY* 

JSH PUSH 


LDY -SXK 

•SXK* 

JSH PUSH 


LOA -FMUL 


JSR C««40 

AY»SXK 

LOY ^2 


JSR POP 

XT 'OZ* 

CALCULATE 

•U3* 

LOY -AZ 

•AZ* 

JSR PUSH 


LOY -BY 

•BY* 

JSR PUSH 


LOA -FRUL 


JSR CMNO 

AZ«BY 

LOY -AY 

•AY* 

JSR PUSH 


LDY ^Z 

•BZ • 

JSR PUSH 


LOA -FMUL 


JSR CMNO 

AY«BZ 

LOA -FSUB 


JS.. Ct«0 

AZ*BY-AY»BZ 

LOY -U3 

LXATION FOR *U3* 
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ORIGINAL PAGE IS 
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JSM HOP 


i 





I 


i 


I 


1 

! 

ri 







CALCULATt 

•UU* 

OKIGINAL PAGE IS 

LD7 -Ul 

•uc 

OF POOR QUALITY 

JSH POSH 
toy 

'UP 

JSR PUSH 
LOA -F»»OL 
JSR Ct«0 

0l*UI 


LOT -02 

•U2' 


JSR PUSH 
LOY ^J2 
JSR PUSH 
(.DA aFMUL 
JSR OIRO 

U2*02 


LOA •*/®0 
JSR OfO 

Ut»UI*U2»U2 


LOT -UO 

LOCATION FOR 'OU' 



JSR POP 


CALCULATE •CDOY* 


LOY -UU 
JSR PUSH 

•UU' 

LOY -03 
JSR PUSH 

•U3' 

LOY -U3 
JSR POSH 
LOA -FHJL 

•U3' 

JSR CHNO 
LOA -FS(fi 

U3«U3 

JSR CMNO 
1 OA -SORT 

UU-U3**03 

JSR CWIO 

S0RT(UU-U3*U3) 

LOY HI? 
JSR PIjSH 
LDA -F(UL 

•U2' 

JSH ONO 

U2»SgHT(0lMl3»U3) 

LOY -03 
JSR PUSH 

•U3' 

LOY ^1 
JSR POSH 
LOA -FUJL 

•01' 

JSR C»MO 
LOA -FAOO 

U3*UI 

JSR C*«« 

U5*U1*U2»S<JHT(UU-U3*U3) 

LOY -UU 
JSR PUSH 
LOA -FPIV 

•OU' 

JSH CW«0 

(U3*Ulm2*SgHT(UU-U3*U3) )/UU 

LOY -CDBY 
JSH POP 

LXATION FOR 'COeY' 

CALCULATt 

•Tf»<S' 

LOY -AY 
JSR PUSH 

•AY' 

LOY -BY 
JSR PUSH 

•BY' 

LOY -COBY 
JSR PUSH 
LOA -FAOO 

•COBY' 

JSR CMIO 
LOA -FDIY 

BY*COeY 

JSR CMNO 
LOA -ATAN 

AY/(8Y*C08YJ 

JSR OMO 

ATAN(AY/(BY«C08Y) » 

LOY -TKMS 

LXATION FOR 'TIMS' 
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a 



I 


i 



M 


97 


JJR POP 


CALCULATE 


LOT 

JSR POSH 
LOA -OOS 
JSR CMNO 
LOT CB 
JSR POP 


ORIGINAL page IS 
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COS(Trt<S> 
LOCATION FOR *C8* 


CALCULATE »CA* 




LOY -THMS 
JSR PUSH 
LOY ^Y 
JSR PUSH 
LOA 

JSR CMNO 
LCM -COS 
JSR OmO 
LOY -CA 
JSR POP 


LOY •TmS 
JSR PUSH 
LOY M>Z 
JSR PUSH 
LOA •hfCO 
JSR C^O 
LOA -COS 
JSR C»t40 
LOY -CC 
JSR POP 


T^S^PY 

COS(T^♦4S^Y> 
LXATION FOR 'CA' 


CALa'LATE 'CC» 


T^f4S-»PZ 

00S(T^^<S4PZ) 
LOCATION FOR 'CC* 


CALCULATE »F» 


LOY -Cl 
JSR PUSH 
LOY -CA 
JSR PUSH 
LOA -FMUL 
JSR CMNO 
LOY -C2 
JSR PUSH 
LOY <0 
JSR PUSH 
LOA -FMUL 
JSR CMNO 
LQA -FAOO 
JSR CMNO 
LOY -C5 
JSR PUSH 
LOY -CC 
JSR PUSH 
LOA -FMUL 
JSR CMNO 
LOA -FAOD 
JSR CMNO 
LOY -F 
JSR POP 


Cl*CA4C2«Cb 

•C3* 


ci*CA4c:=a:40«cc 

LXATION FOP *F» 


calculate »G* 


LOY -C4 
JSR PUSH 
LOY <A 


- 98 - 




i 


I 

i: 


i 


i! 

H 


JSR 

PUSH 


LOA 

•FMUL 


JSR 

CMNO 

C4«CA 

LOT 

■C5 

•C5» 

JSR 

PUSH 


LOY 

■C0 

•CB* 

JSR 

PUSH 


LDA 

•fhul 


JSR 

OMD 

C5«CB 

LOA 

-FADD 


JSR 

CMNO 

C4*»CAK5*C0 

LOY 

<6 

•C6» 

JSR 

PUSH 


LOY 

-cc 

•CC* 

JSR 

PUSH 


LOA 

-FMUL 


JSR 

CWO 

C6*CC 

LOA 

•F^O 


JSR 

CMNO 

C4«CA4C5«C&K:6**CC 

LOY 

■G 

•G* 

JSR 

POP 

GET *G* 


CALCULATE 

•H* 

LOY 

-C7 

•C7* 

JSR 

PUSH 


LOY 

-CA 

•CA* 

JSR 

PUSH 


LOA 

-FMUL 


JSR 

CMNO 

C7**CA 

lOY 


•C8* 

JSR 

PUSH 


LOY 

-CB 

•CB* 

JSR 

PUSH 


LOA 

-FItJL 


JSR 

CMKO 

C8**CB 

LDA 

-FAOD 


JSR 

cmo 

C7**CA4C8**C8 

LOY 

•Ci 

•C9* 

JSR 

PUSH 


LOY 

-CC 

•CC* 

JSR 

PUSH 


LOA 

-FMUL 


JSR 

CMNO 

C9*»CC 

LOA 

-FADO 


JSR 

cm) 

C7**CA4C0**C04C9«OC 

LOY 


LOCATION *H* 

JSR 

POP 



CALCULATE 

•THGS* 

LOY 

k; 

•G* 

JSR 

PUSH 


LOY 

-CIO 

•CIO* 

JSR 

PUSH 


LDA 

-FADO 


JSR 

CMNO 

G4C10 

LOY 

-F 

•F* 

JSR 

PUSH 


LOY 

-cn 

•cn* 

JSR 

PUSH 


LOA 

-f/«0 


JSR 

CMNO 

F4CI 1 

LDA 

-FOIV 


JSR 

CMNO 

(G4CIO)/(F>Cn) 

LOA 

-A TAN 


JSR 

CMNO 

ATAN((G^C10)/(FK 

LOA 

-PTOf 


JSR 

CMNO 

DUPLICATE STACK 


ORIGINAL PAGE IS 
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1 

1 


1 


] 
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L 




t 





« 




LOT "THCS 
JSR POP 

«.T THCS 

LOY -<>tdO 

180/PI 

JSR PUSH 
LDA ■FHJL 
JSR CMNO 

CONVERT FROI RADIANS TO DEGREES 

LOK -3 

POINT TO LONGITUDE FIELD 

LOA -6 

IIOEX LIMIT 

STA XLIM 
LOA -S80 

VIDEO LOCATION FOR LONG. 

STA ¥Y 
JSR T0BCD2 

CONVERT TO DEGREES, MINUTES, SECONDS 

CALCULATE 

•PHGS' 

LOY -THCS 

•TH6S' 

JSR PUSH 
LOA -SIN 
JSR CMNO 

SIN(THGS) 

LOY -CM 
JSR PUSH 
LOY -Cl 4 
JSR PUSH 
LOA -FMIL 
JSR CMNO 

CU«CI4 

'.OA -FMUL 
JSR CI44D 

CI4»CI4»SIN(THGS) 

LOY -H 

•H' 

JbR PUSH 
LOY <12 

•CI2' 

JSR PUSH 
t:‘ FWD 
JSR CMNO 

mci2 

LOA -FMUL 
JSR OMO 

CI4«CI4»SINITICS)*(H+I2) 

LOY < 

•G' 

JSR PUSH 
LOY <10 

•CIO' 

JSR PUSH 
LOA -FM)0 
JSR CMNO 

0<I3 

LOA -FOIV 
JSR C*MO 

CI4»CI4»SIN(THGS)»(H<12)/UXI0) 

LliA -ATAN 
JSR C*MO 

ARCTAN(CI4»CI4*SIN(THGS)*(H<I2)/JG<1, 

LDA -PTOF 
JSR CMNO 

DUPLICATE STACK LXATIONS 

LOY -PHGS 
JSR POP 

GET PHGS 

LOY -P180 

leo/Pi 

JSR PUSH 
LDA -FHJL 
JSR CMNO 

CONVERT PHGS 1HH4 RADIANS TO DEGREES 

LOX -0 

LXATION FOR LATITUDE FIELD 

LDA -5 

INDEX LIMIT 

STA XLIM 
tOA -J60 

VIDEO LXATION FOR LAT. 

STA VY 
JSR TQBGD2 

CONVERT TO DEGREES, MINUTES, SECOICS 


ORIGINAL PAGE jS 
OF POOR QUALITY 


RMO-THETA COWERSION 

calculate *CFI* 


■I 


LOT -PHGS 
JSR PUSH 
lOY -PHI 
JSR PUSH 
LOA -FAOO 
JSR Cl«40 


•PHGS* 

»PMI» 


PHGS^HI 
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t 






t 


LOT -TWO 
JSH PUSH 
LOA -roiv 
JSR cm> 
LOA <OS 
JSR cmo 
LOT -CM 
JSR POP 


r" 


t$2 •• 


(PMGStPMI)/2 



ORIGINAL PAGE !S 
OF POOR QUALITY 


i 

■( 


COS( (PMCS^PMU/J) 
LKATION for 'CFP 


CALCU.ATI 'OLAM 


LOT -THCS 
JSR PUSH 

LOr -RL/W* 

JSR PUSH 
LOA -rsuB 
JSR CW40 
LOY OLAH 
JSR "OP 


•THCS* 

•RLAR* 


THGS-RLAM 

LOCATION FOR •DLAN' 
GCT IT 


CALCULATE '0(411' 


LOY -PHI 
JSR PUSH 
LOY H»HGS 
JSR PUSH 
LOA -FSUB 
JSR CIMO 
LOY -DPHI 
JSR POP 


'PHI' 

'PHIS' 


l -PHCS 

I.XATION FOR 'OPHI' 


I 


CALCULATE 'ARC' 


LOY 

^LAN 

'OLAM' 

JSR 

PUSH 


LOY 

-CFI 

'CFI' 

JSR 

PUSH 


LOA 

-FHJL 


JSR 


DLWPCF 1 

LOA 

^OF 


JSR 

CMNO 

DUPLICATE STACK LOCATIONS 

LOA 

»fnul 


JSR 

Cl«(0 

(0LAM»CFI)«»2 

LOY 

cO«H 

'DPMI' 

JSR 

PUSH 


LOY 

•DPMI 


JSR 

PUSH 

'UPHI' 

LOA 

■FUJI 


JSR 

C(««0 

0PM|*0(4<I 

lOA 

-FJOO 


JSR 

CMNO 

(OLAM»CFI» •'7«0PMI »0HHI 

LOA 

•SORT 


JSR 

CMK) 

SOHT( ... ) 

LOY 

•RCR 

'RW 

JSH 

PUSH 


LOA 

•FMUL 


JSR 

CMNO 

RCR*SgRTl .„. ) 

LOX 

•0 

POINT TO RANGE FIELD 

LOA 

•»eo 

VIOeO LKATION FUR RANGE 

STA 

VY 


JSR 

RNQB 

GE T RANGE 1 N BCD 


CALCULATE 

'PSI ' 

LOY 

■DLAM 

'OLAM' 

JSR 

PUSH 


LOY 

<FI 

'CFI' 

JSR 

PUSH 


LOA 

•f(t(L 



4 



•ll 
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L30 


i 

l 33 






JSR 

cm) 

OtAM»CFI 

LOT 

•OAHI 

•OPHI' 

JSR 

PUSH 


LOA 

•PDIV 


JSR 

CMNO 

OLAN»(yi/OPHI 

LOA 

■ATAN 


JSR 

CMNO 

ATAN(OL^CFI/DPH 1 ) 


ADJUST »PSI 

• FOR mOPER QUADRANT 

LOT 

H)PHI 


LOA 

(BASE),Y 

GET TIC SIGN OF *0»AII' 

BMI 

L30 

GO IF NEGATIVE 

LOV 

-OLAM 


LHA 

(BAS£),Y 

GET SIGN OF 'DLAH* 

0RL 

L33 

GO IF POSITIVE 

LOA 

^UPI 


JSR 

CI940 

PUSH PI 

cDA 

^Pl 


JSR 

CMNO 

PUSH PI 

ItA 

-FAOO 


JSR 

cm) 

ioltiply PI by TNO 

LOA 

•ffCO 


JSR 

cm) 

PSI ♦ 2«PI 

JMP 

L33 


LOA 

•^UPI 


JSR 

CMNO 

PL5H PI 

LOA 

-FADC 


JSR 

CMNO 

PSIAPI 

LOT 

^180 

leo/Pi 

JSR 

PUSH 


LOA 

-FitJL 


JSR 

cmo 

CCNVERT TO OEWEES 

LOX 

-2 

LOCATION FOR BEARING 

LOA 

-SOO 

VIDEO LOCATION FOR BEARING 

STA 

Vt 


JSR 

RNGB 



DONE 



seu SET DECIMAL MODE FOR SENSOR 

RTS 


OWG<r4Al. 

OF POOR QUAL' 


this SUBROUTINl INITIALIZeS T»C PIA FOK USt MITH 
TMi: HATH CHIP AND ThFN SETS THE OONTHOL INPUTS OF T>€ 
HATH CHIP TO INACTIVE STATES. 


.i>A 

PIAA 

1 OA 

PIAB 

IDA 

•$14 

:iA 

"0 

> 

> 

♦ 

LOA 

•0 

sta 

PI AB^t 

LOA 

•SCF 

STA 

PIAB 

LOA 

•4 

sta 

PIAB^I 

LOA 

•7 

sta 

AOCB 

STA 

PIAB 

LOA 

-SF 


CLEAR INTEHHUPTS 

SET INTEKRUPT CONTROL AND DOR 
SET OCR LON 

SET l»#MTS AND OUTPUTS FOR 9115 
SET dor bit HISH 

SET BAOGROUNU coordination BYTE 
SET CO. RD A«l NR HIGH 
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MVE1 


l*USH 


PSH1 


STA P|AB 
LDA PIAA 
LOA -0 
STA BASE 
LOA -2 
STA OASEtl 
LOA -lA 
STA OVSR 


SET SVACK HIGH 
CLEAR ANY INTERRUPTS 
GET LSn OF NIXBER TABLE 
SAVE FOR INDEXING 
CE^ MSB 


ORIGINAL PAGE IG 
OF POOR QUALITY 


STORE 10 FOR BINARv DIVISION 
MOVE NUMBER TABLE TO RE AD/NR I TE SPACE 


LOY 

-0 

CLEAR 

INDEX ONE 

LOX 

■0 

CLEAR 

IK)EX TWO 

LOA 

TABLE 



STA 

BASE1 



LOA 

TABLED 1 



ST'. 

BASE HI 



LOA 

(BASE1 ),Y 

GET A 

DIGIT 

STA 

S200,X 

STORE 

IT IN NEW LOCATION 

ir« 




INY 




CPY 

•120 

MOVE 

120 YETV 

BMI 

MVE1 

\f NOT, CONTINUE 


INITIALIZE VIDEO DISPLAY OUTPUT 


LOX -0 

WRSCN LDY SCLC,X 
LOA lSCRN,X 
AND -sy 
STA S2iOO,Y 
ir« 

CPX -16 
BMI WRSCN 

RTS 


S.BROinriNE TO SEIO DATA TO 9511. 

IT.74 J05BER IS IN REG-Y. BASE ADDRESS 
IS '3ASE.» 


JSR 

SAO 

INY 


INf 


INY 


LOX 

•4 

LOA 

(BASE) ,Y 

STA 

PIV. 

LDA 

•SA 

JSR 

SVIO 

LOA 

PIAA 

DEY 


OEX 


0NE 

PSM1 

RTS 



SET PIAA TO OUTPUTS 
ADJUST REG-Y SO IT 
POINTS TO LS8 OF NUMBER 

LO/J COUNT OF 4 

GET A BYTE OF TRE NUMBER 

GIVE IT TO 951 1 


CLE/V^ ANY INTERRUPTS 
NEXT BYTE 

LOOP UNTIL ENTIRE WORD WRITTEN 
IF DONE, return 


SUB'wUTINE SETS UP PI A SIDE A AS OUTPUTS. 



SAO 


PIAAM 


ca Tl€ CONTROL REGfSTER 
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r 

r • 

f 

n 

f fMC -tFB SET ACCESS THE OCR MAI Df. -ir IS 

1 STAPIAA^I RETURN IT ORIC... »AL K/. « 

i LOA -IFF OF POOR QUALITY 

W %T^ PIAA SET ALL TO OUTPUTS 

^ ^ * 

' * 

K LOA PIAA^I 

ORA -4 SET DOR BIT HIGH 

E STA PIAA^I RETLFN IT 

H 


If?. • THIS SUBROUTINE POPS A NiWbER OfF OF TH£ • 

It • 9511 STACK. NUMBER IS RETURKED TO LXpsTION • 

M • N»TM 'BASE* AS BASE ADDRESS; ITEM IN Y. • 

^K- 



K POP JSR SAI SET PIA AS IK>UTS 

K. LOX -4 LOAD COUNT OF 4 

K POPI LOA >9 


m JSR SV9 

IT STA (BASE),Y STOHE IN TABLE 

|| LOA -SB 

m STA AOCB 


E ORA AOCF 

M STA r'lAB SET RU HICH TO liCK. STACX POINTEH 


K 

K BNE POPI DO 4 BYTES 

H 

W RTS 



*v^H 



W • SUBRr.UTlNE SETS UP PIAA AS INPUTS. • 




fm SAI LOA PIAA<H 

K A^0 -SFe SET OCR BIT LOW 

■k STA PI AMI 

n 

^ LOA -0 

B STA PiAA SET SIDE A TO ALL l»f*UTS 


Wl LDA PIAMI 

V ORA M SET DIF BIT HIGH 

K STA PIAA41 

^ n 

P : 

R'.T' 

B ■ SUBRDUTIFC sends the command byte IN ACCUH. * 

K • TO THE 9511 FOR EXECUTION. ROUTINE PEUFNS TO • 

• caller REGARDLESS OF WHETHER EXECUTION IS COMPLETED • 

Hf • OR NOT. • 

I- 


^ te 

F^A 

H JSR SAO SET PIA SIDE A AS OUTPUTS 

H pla qet tm command 

R STA PIAA SEND TO 9511 

K LOA -IE 

b JSR SVIO 


BIT PI/'JL*! TEST IF COMMAK) DONE 


HSE 8PL •-> KEEP testing until DOt« 

IK LOA PIAA CLEAR THE INTERRUPT BIT 

s 

|E • REM) THE STATUS REGISTER; RETlAtH IF INVALIU CODE 

9 

R 


E 


ft 

fli 


■iv - 

B -104- 

m-. 

1 
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MAS PROOUCEO. 


ORIGINAL PAG^ i:, 
OF POOR QUALiTY 


i 


I 



I 







i 


r 


t 




(X 


JSR SAI 
l0^ -SO 
JSR SV9 
PHA 

LOA -SF 
STA AOCB 
ORA AOCF 
STA PIAB 
PLA 

AND -noon no 

BEO OK 
PLA 
PLA 
SFI) 

RTS 


^^ET PIAA SIDE A AS INPUT*: 

SET C/D HIGH, RD LON 

SAVE IT 

RETURN 9511 TO INACTIVE STATE 


ZERO OUT UNIMPORTANT bITS 
IF ZERO. CONTINUE PROCESSING 
POP THE STACK 

SET DECIMAL FOR SENSOR ROUTINE 
RETURN 


STA 

AOCB 

ORA 

AGCF 

STA 

PIAB 

BIT 

PIAA-fl 

BVC 

•-J 

LOA 

PIAA 

RTS 



SVIO 


PHA 

STA AOCB 
ORA AGCF 
STA PIAB 
PLA 
CLC 
ADC -1 
STA AGC0 
ORA AGCF 
STA PIAB 
RTS 


miS routine CONVERTS THE BINARY FLOATING-POINT LATITUDE 
AND LONGITUDE SEPARATELY INTO THE STANDARD DEGREE. MIN.ITES, 
AND SECOK) FORMAT. THE RESULT IS STORED IN BCD. 


INTG CLD 

^TOF 
JSR CMNO 
LGA -FIXO 
JSR CMNO 
LOA -PTOF 
JSR CMNO 
loy -temp 
STX XTEMP 
JSR POP 
DEY 

LDA (BASE).Y 
STA 0V0N4-1 
DEY 

w'X (BAS£),Y 
STA OVON 


SET decimal MODE OFF 

DUPLICATE ST^CK LOCATIONS 

CONVERT POSITION TO AN INTEGER 

DUPLICATE IT 
EORK AREA 


GET T^€ HEX RESULT 

AND PLACE FOR HEX-TO-BCD CONVERSION 


ROUTINE TO CONVERT A TWO-BYTE HEX NIA40ER TO A FIVE-BYTE BCD 
NUMBEF^. 



\ 






1 




■ 

LOX 

•4 

COUNT OF FOUR 

• 

t 


OiYIUe BY 

TtN 

UNSPO 

LOA 

•0 



STA 

Rf«40P 

CLEAR RtNAINOER 


LOY 

-17 

SET UP COUNT 


JNP 

001 


007 

LOA 

RMNOR 



sec 


SET CARRY FOR SUBTRACT 


SBC 

OVSR 



BPL 

NREST 

00 IF NO RESTORE 

001 

CLC 




JNP 

MERGQ 

GO TO SET g 

N«EST 

STA 

RI#IOR 

NEW HESIOUE 


sec 



MCRGQ 

ROl 

OVO»H1 



ROC 

OVON 



OEY 


OECKENENT CXXJNI 


eeg 

RTN 



KOI 

RN^iOR 

SHIFT LEFT 

• 

JNP 

002 

CONTINUE 

KTN 

LOA 

RNNOR 

GET HEmiNOER 


STA 

ONES^X 

STORE IT 


OCX 



• 

BPL 

UNSPO 

00 UNTIL DONE 


LOX 

XT£>f> 


• 

t 

RTS 



T0eCD7 

JSR 

INTG 

GET M. INTEGER PART 




ORIGINAL PAGE IS 
OF POOR QUALITY 


STOKE THE HOSITION COOKOINATES ON ^lOEO SCKEEN 


LOY 

VY 

GET VIDEO LOCATION 

LOA 

ORESO 

LOM) MSII 

ORA 

-ISO 

CHANCE TO ASCI 1 

STA 

S2200,Y 


INY 

LOA 

ORES-M 

Of LSfl 

ORA 

-$50 


STA 

12200. Y 


INY 

INY 

STY 

VY 


LOA 

0RES43 


ASL 

A 

MOVE LONER FOUR BITS TO UPPER FOUR BITS 

ASL 

A 


ASL 

A 


ASL 

A 


STA 

LAT.X 

STORE THE DIGIT IN POSITION FIELD 

LOA 

ORESH 

GET next digit 

ORA 

LAT,X 

MERGE WITH UPPER DIGIT 

STA 

LAT,X 

REPLACE 

\f^ 

CPX 

XL IN 

SEE IF DONE 

HEQ 

0UT2 

return if done 

LOA 

•fltd 


JSR 

CNNO 

C.4ANCE INTEGER TO FLOATING-POINT 

LOA 

•FSUB 


JSR 

CNNO 

SUBTRACT OFF INTEGER PART 

LOY 

<60 

GET CONSTANT "60* 

STX 

XTENP 


JSR 

PUSH 


LOX 

XT»r 


LOA 

•fmi 





:SK cmo 

jNM roao); 

KfS 


MULT in Y KtSllHJl IIY oO 


ORIGINAL PAGE IS 
OF POOR QUALITY 


houtini ooNvtHrs manui 

AM) STCKtS Tnlm. 


AM) UlAMINl) M4UM HI MVK) TO |ia> 


MMV J$M INTO 


UtT IM INTlliiH MAHl 


»ll)tO OUTPUT 


Ll)Y 

VY 


Ll>A 

OKIS^: 


OKA 

•ISO 


STA 

y20o,Y 


INY 

LCIA 

OKLStS 


OKA 

•ISO 


sta 

W200.Y 


iNY 

LUA 

OKISM 


ORA 

•ISO 


SYA 

W200,Y 


INY 

INY 

^V( 

VY 


lOA 

IWS^2 

XU OKIHK UlUl f 

ASl 

A 

MUVt LWIK UlUlT 10 un*tK 

ASL 

A 


ASL 

A 


ASi 

A 


STA 

KMO.X 

PUT IT IN KLSULT 

U>A 

wests 

2NU OKIHK DIGI 1 


KHO,X 

MlKl'4 with last UIUII 

STA 

KHO,X 


INX 

lOA 

wesM 

f IKST UIGIT 

ASL 

A 


ASL 

A 


ASL 

A 


ASl 

A 


STA 

KHO.X 

STONl IT 

lOA 

•FlTl) 


JSK 

L>t40 

GHANU TO FLOAT INLW*0 INI 

ll^ 

•FSU> 


JSK 

CMNO 

SUITTKACT OFF INILGIK 

LOY 

•C2M> 

-2S6** 

STX 

XTU#» 


JW 

KUSH 


U>A 

•^niL 


JSK 

IXNO 

GtT MIACT K>N AS INTIGIK 

lOA 

•f 1X0 


JSK 

C 


lOY 

•TtMP 

MA»#L AKIA 

JSK 

POP 


OtY 

LOA 

UlASn.Y 

GLT THl FKACIION 

LSK 

A 

SFAKl>4-TAHLt Ml THiIl OF 

lsk 

A 

LUTAININD COKK(SK)Nl)l NL. 

LSK 

A 

UlCiMAL VAIUI 

ISK 

A 


TAX 


list AS 1 NIX X 

LOA 

FKTBLI ,X 

OiT TNI OICIMAl tgUIVALlNl 


- 10 /- 


■Mm 


LOX XTCM> 

LOT VY 
ORA -»0 
STA 12200, Y 
AfC -IF 

ORA RHO.X STORE T»C FRACTION 

STA RHO.X 

RTS 


ORIGINAL PAGE IS 
OF POOR QUALITY 


CONSTANT TABLE OF NIMBERS USED IN CALCULATIONS. 


▼ABLE 

CUTS 


t 



(CNTS) 




EQU 

• 




HEX 

7C,9f .BE.77 

00 

- 

TCY 

HEX 

7C,00,2f ,1B 

04 

- 

TCZ 

HEX 

7E,9A,EC,71 

08 

- 

T^61Y 

HEX 

7E,96,0F,39 

OC 

- 

jmz 

HD 

06,B8,F2,6E 

10 

- 

XNR 

HD 

OO.FD, 13,63 

14 

- 

CTMY 

HD 

7E,9A,55,50 

18 

- 

S'WY 

HD 

00,FD,2E,C3 

IC 

- 

CIHZ 

HEX 

7E, 97, 60, 51 

20 

- 

STMZ 

HD 

7F,BF,D3,EB 

24 

- 

CXK 

HD 

00,E0,5A,69 

26 

- 

SXK 

HD 

02, BE, EC, OD 

2C 

- 

Cl 

HD 

03,69,F5,17 

30 

- 

C2 

HD 

83,E5,7E,A9 

34 

- 

C3 

HD 

03,62, 0E,4C 

30 

- 

C4 

HD 

03,6E,OC,13 

X 

- 

C5 

HEX 

01,93,27,AA 

40 

- 

C6 

HD 

03,A1,A2,61 

44 

- 

C7 

HD 

03,A3,CA,CE 

46 

- 

C0 

HD 

00,6F,26,B5 

4C 

- 

C9 

HD 

77,CF,C0,08 

50 - 

CIO 

HD 

75,67,75,21 

54 - 

Cl 1 

HD 

F6,C2,1E,A6 

50 - 

C12 

HD 

01, 60, 6F, 75 

5C - 

C14 

HD 

69,06, BF, 95 

60 - 

IE-7 

HD 

X,06,F5,7B 

64 

- 

RCR 

HEX 

02,80,00,00 

60 

- 

TWO ("2") 

HD 

09,60,00,00 

6C 

- 

C256 ("256-) 

HD 

06,E5,2E,E1 

70 

- 

P10O (160/PI) 

HD 

06.F0, 00,00 

74 

- 

C60 (60) 


• 6N0 Ot NUMUER TABLE; START SCRATCH S»>ACE 


KRTbU HEX 00,01.01.0^,05,05,04,04 
HtX 05,06,06,07.00,06,09,09 
LSCRN ASC •LAlLONCRNQBRfC. .» 
sac HEX 66,69, 6A 

HEX 66,69,8A,8B 
HEX A0,A9,AA 
HEX C6,C9,CA,C8 
HEX B5.03 


END 


108 


original 


Appeadix IV . Program listing for program to com- 
pute spherical C-constant from equations 3-22 through 
3-36. This program was written in FORTRAN-IV using double 
precision arithmetic. 


- 109 - 


o’^no oooo ooo ooo oooooooooo 




this pwograh is used to calculate the constants used in the 
EXPLICIT coordinate CONVERTER TO CONVERT THE GEOCENTRIC 
ARC ANGLES TO GEOCENTRIC JOOROt NATES (LAT. AND LONG.) 

REFERENCE H. FELL, NAVIGATION, SUMMER 1975 AND PRIVATE 
COMMUNICATION. 

FIRST WRITTEN I 1/1980 
MAJOR REVISION 6/1961 
J. P. FISCHER 

IMPLICITREAL»8(A-H,0-Z) 

DATA RR/6578.2064/,R/6356.5838/ 

CCMMON XX,vy,ZZ 
PRINT23 

CALL RDLL$(PHI0,THET0,49) 

CALCULATION OF THE CENTER Of * OSCULAriNG SPHERE. 

BI-OSIN(PHIO) 

B2"BI»B1 
B3-DC0S(PHI0) 

B4-83»B3 
B5-RR»RR 
B6-R*R 

RC"<B5*B2Yfl6«84)/R 
B7»OS0RT(B5*B2*B6»B4) 

XX-(RR-e7)«B3*DSIN(THET0) 

TV- ( RR-B7 ) »DC0S ( THETO) «B3 
ZZ-(R-RR/R»B7)»0I 
CIO-XX/RC 
Cl l■YY/RC 
CI2-ZZ/RC 
CI>CIO 
CU-RR/R 

READ IN THE COORDlfWTES OF THE Th«EE LORAN-C STATIONS. 

CALL R0LLS(PM,TM,A9) 

CALL R0LLJ(PX,TX,49) 

CALL RDLL$(PY,TY,49) 

CONVERT THE GEODETIC LATITUDE TO GEOCENTRIC AND MAP THE POINTS 
ONTO THE OSCULATING SPHERE. 

CALL SPO$(PM,TM,PV*<,Tffl) 

CALL SPO»(PX,TX,PHX,T^«) 

CALL SPOtlPY,TY,PHY,THY) 


ORIGINAL PAGE 13 
OF POOR QUALITY 


COMPUTE THE C-OONSTANTS USED TO CONVERT GEOCEN^IC ARC ANGLES 
TO GEOCENTRIC LATITUDE AND LONGITUDE. 


C 


AI-XOSCP^«<) 

A2-0SIN(PHM) 

A3-0C0S(T^W) 

A4-OSIN(Tl«4) 

A5-0C0S(PHX) 

A6-0SIN(PHX) 

A7-DC0S(THX) 

A8^SIN(THX) 

A9-0C0S(PHY) 

AKW)SIN(PHY) 

A1 l-OCOS(THY) 

AI2»0SIN(THY) 


0-A6»Al«A9*(A3*AI2-A4»Al l)-A2*A5*A9»(A7*AI2-A6*AI I )+AI0*A5* A1 • (A7« 
>A4-A8*A3) 

CI-(AI«A4*AI0-A2*A9»AI2)/D 
CX A6*A9» At 2-A5»A8»AI0)/D 
C > ( A5»A8»A2-A6* Al • A4 ) /D 
C4-(A2*A9*AI l-AI»A3*AI0)/D 
C5-(A5*A7»AIO-A6»A9«AI D/D 
C6«( A6*A1 •A3-A5*A7« A2 )/U 




-no 


ooo oonooo non 


C7«A>»A9»(AJ*Ai;-M»AI l)/0 

Ce»A5*A9»(A8«An-A7«AI2)/D ^ 

C9-A5«AI»(A7«A4-Ae»A5)/0 ORIGINAL PAGc 13 

PRINT THE RESULTS AM) QUIT. OF POOR QUAHTY 

PRINTI,CI,C3,CJ,C4,C5,Ob,C7,Ce,C»,C10,CI l,CI 2 ,C) 3 ,CU 
9 STOP 

I FORMATC Cl - ' ,IPOM.7,3X,*C2 • ' ,014.7, 3X, 'C3 ■ ',014.7/' C4 • ' 
>,0U.7,3X,'C5 • ' ,014.7, 5X,'C6 - ',014.7/' C7 - ' ,014.7, 3X, 'C8 - ' 
>,0I4.7,3X,'C9 - ',014.7/' CIO - ' ,014. 7, 5X, 'Cl 1 - ' ,014.7, 3X,'C12 
>- '.014.7/' C13 • ' ,014.7, 3X,'C14 - ',014.7) 

23 FORMATC '1 ENTER COOROIMLTES OF MIOPOINT Of COVERAGE ARI>. ' > 

END 

SUBROUTINE SPOSCPHIO.THEO.PH 1 1 ,THEI ) 

THIS subroutine calculates the COOROINATES of a r-OINT ON 
A SPHEROID ANP CALCULATES THE IMAGE OF I* UNDER THE SPHEROID 
TO the osculating sphere KITH CENTER AT XX. VY, Zl FROM 
MAIN PROORAM. 

implicit REAL'aCA-M.O-Z) 

DATA RIV6578.2064/,R/6 356.5a38/ 

COMMON XX,TY,Z2 
PHIO-OATAN(R/RR*DTAN(PHIO) ) 

X-RR«OOOS(PH IO)*OSIN(THEO) 

Y-RR*DOOS(PHIO)''KOS(THEO) 

ZM««OSIN(PHIO) 

XL-OS0ftTI(X-XX)»(X-XX)4(Y-YY)»(Y-YY)4(2-ZZ)»(2-ZZ)) 

PMII-DARSIN((Z-ZZ)/XL) 

THEI-OATANC (X-XX)/( Y-Y Y) ) 

REnUN 

END 

THIS SUBROUT I )€ CONVERTS GEXENTP.C COTROINATES ENTERED BY THE 
USER TO RM)IAN COOROINATES. IM>UT FORM IS: 0000 MM SS.SS 
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Appendix V . Prograoi listings for service routine 
to convert BCD numbers entered at user's terminal to 9511 
equivalent binary floating-point numbers. First p. ogram 
is written in FORTRAN-IV and second program is written in 
IBH 360/370 assembly language. 



ooo ooo oooo 


WHEHE '0000' IS THE OECHEES PORTION Of THE LAT. OH LONG., INCLUOINC 
SIGN. MM IS THE MINUTES PORTION. AND SS.SS IS IHE SECONDS PORTION. 
READ FORMAT IS: 14 . 1 X, 12. 1 X.FS.O. 

IMPLICIT REAL*0(A-H,O-Z) 

DATA Pl/3.141?926555090/ 

DATA MSGI/'LATI'/,MSG2/'TUDE'/.MSG5/' : '/,MS04/' LONG'/.MSCV I TU 

>0'/,MSG6/'E: '/ 

PII-PI/100. 

PROMPT USER. 

PR I NT I, MSG I, MSG2, MSG) 

READ(5,IO,ENO-2I)TI ,T2,T3 
PHI-PH«(TU(T2* T3/60. )/60.) 

PROMPT USER FOR LONGI TUOE ENTRY. 

PR I •fl’ t , MSG4 , MSGS , I4SG6 
READ(5.10,END»2I)T»,T2,T3 
»H£*P 1 1 •(TU(T2^T3/60. )/60. ) 

return 

21 RETURN \ 

I FORMATC ENTEI'. '.3A4/' 0000 MM SS.SS') 

10 FORMAT(F4.0,)X,F2.0, 1X.F5.0) 

END 
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C FORTRAN PROGRAM 9)11 

C VHIS PROGRAM IS DESIGNED TO RUN IN CONJUCTION WITH SUBROUTINE 

C DECHX. A REAL NIMBER IS READ IN AND CONVERTED FR(M 360/370 

C FLOATING-POINT FORMAT TO THE BINARY FLOATING-POINT FORMAT 

C USED BY THE AM931IA MATH CHIP. TIC INPUT NIMBER IS HANDLED 

C AS REAL'S TO ALLOW ROUNDING UP OF THE FII«^L RESULT. 

C 

C J.P. FISCHER 4/81 
C 

REAI'0 XA 

LOGICAL'! R1,R2,R),R4 

1 REAO(3,!0,EN>2)XA 
CALL DECHX ( XA ,R 1 ,R2,R3,R4) 

PRINT I I,R1,R2,R3,R4 
GOTO! 

2 STOP 

10 F0RMATCFI5.0) 

11 F0RMAT(IX,Z2.IX,3Z2) 

END 


THIS FOR-mAH-CX)HPATIBLE SUBROUTINE TRANSLATES A 360/370 
FLOATING-POINT NLMBER TO THE FLOATING-POINT FORMAT USED 
BY THE AH9511 MATH CHIP. THIS IS DESIGNED AS / DEVELOPMENT 
PRC3GRAM. RESULTING H\HBER IS RETURNED AS FOUR SEPARATE 
NIA4BERS (LOGICAL^I). 

J. P. FISCHER 3/81 
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OED« 


POSM 


LEFT 


NOT 


POS 


SPACE 

2 


CSECT 

BC 


BRANCH AROUND ID 

DC 

XLP5* 


DC 

CL5* DECHX* 


LR 

12 J3 

NEW BASE 

USING DECHX, 12 


STM 

14,12,12(13) 

SAVE REGISTERS 

SPACE 

MVI 

FLAG,0 

CLEAR FLAG 

L 

2,0(1) 

ADDRESS OF FLOATING-POINT NUMBER 

MYC 

EXPO(5),0(2) 

MOVE INTO THIS PROGRAM 

TM 

EXP0,X'80* 

SEE IF MANTISSA IS NEGATIVE 

B2 

POSM 

00 IF POSITIVE 

MVI 

flag,i 

SET MINUS BIT 

Nl 

EXP0,X*7F» 

CLEAR NEGATIVE BIT 

SLR 

3,3 

ZERO A REGiSTEH 

1C 

3, EXPO 

GET THE EXPONENT 

SL 

3,F64 

SUBTRACT X»40* FOR NEW Blf^RY EXPONENT 

SLL 

3,2 

CONVERT EXPONENT FROM HEX TO 01 

SPACE 

TM 

EXP0^1,X*80* 

SEE IF left adjustment NECESSARY 

BO 

NOT 

IF SET, DON'T DO ANY MORE 

ICM 

4,15,EXPOf1 

GET THE MANTISSA 

SLL 

4,1 

SHIFT OVER ONE 

BCTR 

3,0 

DECREMENT EXPONENT 

STCH 

4,15,EXP0^1 


D 

LEFT 

SEE IF SHIFTING COMPLETE 

TM 

EXPOM,X*80* 

SEE IF THIS BIT IS SET 

BZ 

'\)S 

IF NOT, DON'T ROUND UP 

ICM 

4,B*0111*,EXP0^1 GET THE MANTISSA 

LA 

4,1(4) 

INCREMENT IT 

STO 

4,B*0111»,EXPOf1 REPLACE IT 

SPACE 

STC 

3, EXPO 

STORE THE NEW EXPONENT 

TM 

flag,i 

SEE IF NEGATIVE 

BZ 

P0S1 


01 

EXPO,X*80* 

SET THE NEGATIVE BIT 

B 

P0S2 
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POS^ 

Nl 

SPACE 

EXP0,X»7F» 

CLEAR THE NEGATIVE BIT 

• 

SPACE 

return the new 

Float iNG-POl NT io^ber 

P0S2 

IN 

3«6,4(1) 

ADDRESSES OF RETURN LOCATIONS 


MVC 

on. 3). EXPO 

MOVE THE NEW EXPONENT 


MVC 

0(1 .4),EXP0^1 

MOVE THE MSO OF MV^TISSA 


MVC 

0(1 .3) .EWHH7 

MOVE THE Ml DOLE OF M# NT ISSA 


MVC 

SPACE 

0(1 ,6) .EXROO 

MOVE THE LSB OF MANTI SA 

• 


RETURN TO CALLER 


SPACE 

LM 

2,12.28(1.; 

RESTORE REGISTERS AS THEY WERE 


MYI 

12(13) ,X*FF» 

INDICATE CX)NTHOL RETURNED TO CALLER 


SR 

13,13 

ZERO return CODE 


0R 

14 

RETURN 


SPACE 

3 



OS 

OF 


F64 

OC 

F*64* 


EXPO 

OS 

XL3 


FLAG 

OS 

XL1 



END 

OECW 
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